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ABSTRACT 


Industry often demands the generation of preprogrammed motion. In fixed 
automation, it becomes necessary many a times to have a mechanism which 
generates a predefined motion with high accuracy and repeatability, and at a high 
speed. Computer controlled robots become unsuitable in such cases owing to their 
relatively lesser motion resolution, low repeatability, and high cost . A 
multi— degree— of —freedom mechanism, composed of cam modulated linkages, is 
proposed as a means | of achieving the above requirements. The simplicity in the 
design and manufacture, and capability to generate precisely controlled motion 
including the generation of non-differentiable paths, makes these mechanisms very 
useful. An attempt is made in the present work to develop a procedure for the 
design of these mechanisms. 

The present work comprises of three parts. In the first part, a strategy 
for the design of cam modulated linkages capable of generating a given 
trajectory, has been developed. Along with a procedure to synthsise and analyse 
the mechanisms, critical hardware and practical design issues are also discussed. 
Examples supporting the proposed method are presented. 

The second part is an extension of the design of cam modulated linkages, in 
which, a method is proposed for the optimal synthesis of cam follower mechanisms, 
to minimize the cam size in cosideration of operational constraints. 

In the third part, software to assist the design of the cam modulated 
linkages incorporating the above two parts has been developed. 



chapter 1 . INTRODUCTION 


IjL Introduction. 

In fixed automation, it is frequently necessary to generate programmed 
motion which remains mchanged for a large number of operation cycles. Typical 
applications include mechanical part handlers on assembly lines, and profile 
cutting machines operating at a mass production rate. The demands on the machines 
producing such a motion are that they should be inexpensive, accufl'ate, and 
preprogrammable. Though computer controlled reprograrranable robots can be used 
in such cases, they have certain limitations in that they are expensive, have lower 
repeatability, and lower motion resolution. Also, th^y cannot be made to operate at 
very high speeds. In such cases, mechanical devices termed as mechanicsl robots, 
can be employed. In these devices, the output motion is governed through 
mechanical elements like 'cams and linkages, and is achieved by pre-designing the 
geometry of these elements. 

Many of these applications generally demand more than one degree-of- 
freedom in the output motion. The profile cutters for example, require a feed rate 
control, artti a spaed or depth of cut rate control, independent of each other. 
Similarly, the assen^ly part handlers may need a stroke and a lift, or stving about 
some axis, and a lift, which again may be independent of each other. There is 
generally one single prime mover in mechanical devices from which power is 
extracted. To have more than a single degree-of -freedom motion, it becomes 
necessary to have as many rKimber of constrained mechanisms as the number of 
degrees-of-freedom, each taking power from the same shaft, and each controlling 



one component of motion. Thus the n dimensional motion gets divided into n 
components, each generated independently by one mechanism according to specified 
requirements. Planar mechanisms can be used to govern the motion in each 
direction, and this proves to be a simpler choice from design and manufacturing 
point of view. 

In as far as the generation of planar motion is corwemed, linkages will be 
the first choice. However, when a precise control over the output position and 
velocity or acceleration throughout the motion is needed, linkages become 
unsuitable. For each specified point of the output motion, there will be one 
mathematical condition in terms of the geometric parameters of the mechanism. 
Therefore, to have an infinite point accurate motion, the mechanism requires an 
infinite number of geometric parameters to specify its dimensions. The cam 
mechanism is such a mechanism; the points on the cam profile being the design 
parameters. However, as is shown in the present work, the use of just a cam' 
mechanism is not sufficient in nrost of the cases, since the range of motion that 
can be extracted out of a cam is generally small, unless the cam itself happens to 
be large. Also, most of the time, it is batter to use an oscillating roller follower 
since it permits the use of a larger pressure angle. In such cases, there wou^ be 
a necessity of converting angular motion to linear. Moreover, in some cases, the 
cam profile may turn out to be too complex to be of any practical use. In all such 
situations, if this higher pair could be used as a part of a cam modulated linkage 
(CML for short), all purposes — viz. magnification of motion, conversion from 
angular to linear (and vice-versa), and alleviating the demands on the cam profile 
through the input output relationships of the linkage - wfxild be served. Hence a 
cam modulated lird<age becomes a better choice in this application. A multi degree- 
of-freedom mechanism with planar cam modulated linkages governing one degree of 



freedom eachj is proposed in the present work. 

A large number of methods for the design of linkages and cams taken as 
individual elements have been developed. <Ref [13, t23, [53, [63, [73, [113) However; only 
a few methods exist for the design of cam modulated linkages. Number synthesis of 
CM_s has been done by Kohli arwl Singh [83, ard Sandor and ErdmeUi [113. Sandor and 
Erdman [113 have discussed the nessecity and sufficiency criterion for a single 
linkage to be corwerted to a cam modulated linkage. A few analytical methods for 
the design of ChLs have been pjroposed by Kohli and Singh [83. 

The fii^t ift¥X5rtant step in the design of CMLs is that of number sijithesis 
in t^ich either by the enumeration and classification techraque [113 or otherwise, 
all CM-s having the same constructional features are evaluated, and a selection 
fr'om these based on the application is made. In our case, it can be seen that the 
cams and linkages used in the mechanism actually serve two different purposes. 
Cams in particular are used to have a programmed output irotion in which the 
displacements, velocities, accelerations and dwells if any, take place in a 
kinematically controlled manner. Linkages on the other hand serve the purpose of 
motion magnification, and conversion of motion from one type to another. (Angular 
to linear, etc 3 This shows that a modularity exists between the linkages and cams. 
Hence it should be easier to design and manufacture these units if cams te used 
in the first stage of motion transfer and linkages in the second. The input from 
the shaft at a constant angular velocity, can be converted to a desired time 
functictfi of the follower displacement, velocity etc. by the cam; and the linkage may 
further magnify or convert this motion to a desired ti^, prior to passing it to 
the output end. Coordination and combination of individual motion components form 
an important part of design arrcl is achieved through elements termed as moiion 


assemblers . 
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Fig. 1.1 shows a schematic sketch of a 2 degrees-of -freedom motion 
generating mechanism based on the above principle^ producing motion in the x-y 
plane as a combination of motions along x and y directions. As can be seen from 
the figurej both the motions are independent of each other, and are finally coupled 
by means of the assembler. 

The second step in the design of mechanisms comprises of dimensional 
synthesis in which geometric parameters of all the components in the system are 
designed. To restrict the size of the cam within reasonable limits, it becomes 
necessary tc limit the cam«follower displacements within a certain range. Linkage 
synthesis primarily comprises of coordinating this range with the end effector 
motion range in that particular direction. 

Cam design comprises of designing the follower parameters, (e.g. the arm 
length, roller radius etc., for an oscillating roller follower) and the synthesis of 
the cam profile. In designing the follower parameters, the approach is to reduce 
the cam size to the maximum possible extent, and an optimizing algorithm has been 
used for this purpose. Cam profile design is based on these follower parameters^ 
and the follower motion which can be determined by performing the inverse 
analysis of the designed linkage. 


12 Scope and abjective of the noHc. 

Industry has been using the type of mechanisms described above. However, a 
formal study regarding the different design issues remained to be done. It was 
sought at the conmencement of this work to explore these aspects, and to inclement 
them in the form of an interactive programme. 
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A generalised approach for cam and linkage desi^ was felt necessary to 
fasciliiaie computation. Most of the methods for linkage design use analytical 
expressions needed for the particular linkage under consideration. Hwever, since 
here, a large number of such linkages would be involved, it would not be feasible 
to use analytical expressions for each of the possible linkages. Also, doing so 
would limit the application of the procedure to only those linkages i^oh were 
considered at the stage of implementation. 

Since the design primarily involves the determination of all the geometric 
parameters of the mechanism, a need was felt to develop a strategy for the design 
of the cam-follow€M' parameters t«^ch are normally based on general judgement, 
considerations of space etc., or deduced from some standard noTMjgrara. 

In view of the above aspects, the following points reflect the objectives of 
the present work. 

— Formulation of the linkage synthesis and analysis problem in a generalised 
form. 

— Development of a design scheme suitable for the mechanisms described 
earlier, by bringing together existing methods of cam and linkage design. 

— A general study as regards the applicability and limitations of U^se 
mechanisms, including some related practical issues. 

— Building tffj a strategy for the optimal desisn of the follower parameters. 

— Developmait of a software incorporating all these ideas, followed by the 
tasting with a few possible mechanisms. 

Scope. 

In the present work only 2 degrees-of-freedom systems have been dealt with, 
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and the resulting motion of two, and two and half dimensions of space (Sect. 4.2.1), 
has been discussed. The principle remains the same for a larger number of 
degrees-of -freedom. The cams considered are of the disc tgpe. The method for 
design however, is developed in such a way that 3-d cams can also be included. 
Linkages are assumed to be of four joints only, with any cwnbination of revolute 
and prismatic joints. In most of the cases these will suffice and it will be rarely 
necessary to go teyond what is offered by a four bar chain for this application. 

Division of the work. 

With the objectives as mentioned above, the work has been subdivided into 
the following parts. 

1. ChC desigri. 

2. Follower parameter optimization. 

3. Software development incorporating (i> and (2) above. 

OIL design has been discussed in chapter s 2, 3, ard 4. The design procedure iu 
developed in Oapter 2. Okapter 3 discusses the details of the computaticwul and 
other details of th% design. Chapter 4 deals with the hardware, and practical 
issues. Chapters 5 and 6 are devoted to the design of the follower parameters in 
view of optimization of the cam size. Details of the implementation appear in 
chapter 7. Finally, the conclusions and possibilities of further work are mentioned 
in chapter 8. 
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chapter 2. DESIGN PROCEDURE 


2J. Preliminaries. 

The design problem under consideration can be staled as: 

"Given an arbitrary Irajeclory expressible in two planar independent 
variables (angular or lirear), to design cam modulated lif^ages which will be 
capable of generating this trajectory." 

The procedure consists of the following steps: 

1. Trajectory specification, in which the path to be followed and the end 
effector kinematics are represented in a way so as to fascilitate the 
SLdDsequent st^thesis. 

2. S^thesis of C^m Modulated Linkages (CHL) for each independent direction, 
ie for each degree of freedan. 

3. Assembly and coordination of the irwDtion generated by the CM- for each 
direction. 

4. Analysis of the designed mechanism to test its performance, and 
resynthesis if needed. 

These are explained in the following para^aphs. 

22 Systems of coordinates. 

We shall in ger«ral have Df * 2 permenatly fixed coordinate systems, (Fig. 

2.i) hdiere Df is the number of CMLs in the mechanism i.e. the number of deff'ees-of- 
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Z-1 System of Coordindes. 


freedom (2 in our case). These will be as follows. 

- One global coordinate system G. 

- One local coordinate system Ej for the expression of end effector 

displacements. 

- D/ fixed local coordinate systems for the description of linkage 

configurations. (One coordinate system per linkage.) 

- Df temporarily fixed coordinate systems These are fixed at the cam 
centre. They are termed as temporary since the location of the cam centre is 
unknown at the commencement of design. The temporary coordinate system is used 
to design the cam-follcsMer mechanism) and after all the cams in the system are 
synthesised, it is transported to a permenant position during the assembly of the 
CMLs. 

In addition to these, we will have local moving coordinate systems used in 
the cam profile si^pithesis; uriiich are explained later. 

23 Trajectory specification. 

The trajectory is the motion generated at the end effector of the 
mechanism. Curves in space are generally represented in terms of D-1 simultaneous 
equations, where D is the number of dimensions of the space trtder consideration. 
In our case, let Xi and Xg be the two dimensions, with time (t) as the third 
dimension. Then the following forms of representations hold. 

1. Explicit specification of X 2 and Xi as functions of t. 

Xi = fi(t), X 2 X f2(t>. 

2. Explicit parametric etjjations. 

Xi = fi(u), X 2 = f 2 (u). t = f 3 (u). 


iO 



3. One space -Function, one space-lime function. 

fitXi, Xa) = 0, faOCi, Xj) = gd). 

Typically, fj will represent the path generated in space, fa the length along' 

the curve, and g, its time variation. 

4. Implicit relations in Xj, Xa and t. 

fiOCi, Xa, t) = 0, faOCi, Xa, l> = 0. 

These are generally not enoountered in practice. 

5. In the form of numerical tables over the cycle time. 

The actual trajectory followed may be a single continuous function, a 
piecewise continuous function, or a discrete numeric function. In any of the above 
forms of representation, the values of Xi, and Xa at any given time instant can be 
determined by the solution of simultaneous equations, or by direct substitution. A 
search in these directions gives the extreme values of the displacement needed 
for lirdcage design. 

It can be seen that since Xi(t), or X 2 <t) are the motions generated by one 
single CML, the higher derivatives of these with respect to time need to be 
continuous for a smooth performance. However, there is no such restriction on 
XiOCa). Thus it is possible to generate motion which can be non-differentiable at 
intermediate points i.e., sharp corners may as well be incorporated in the 
trajectory. 

The trajectory is always represented as a curve in the fixed frame of 
reference E|. 

2.4 Synthesis of 01_s. 

S«,rithesis comprises of the following: 
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1. Linkage synthesis, based on the needed range of the output motion and 
desired range of input motion to the linkage. 

2. Inverse analysis of this linkage, to determine the input motion needed over 
a complete cycle. 

3. Design of the cam follower system based on the input motion to ttw linkage. 

2.4J. Linkage synthesis. 

Objective of the linkage synthesis is to coordinate the required range of 
output motion and desired range of input motion. A generalised nethod for the 
design of four bar linkages is presented. 

Scheme of representation. 

Representation in the form of vectors has been used here. Q^ig. 22) Though 
we have only four link mechanisms under consideration, an n vector model with one 
degree-of -freedom has been used. This is to fascilitate the representation of links 
which can be expressed in the form of more than one vector, eg. a bell crank 
lever. This is closer to the physical solution in the first place, and is especially 
useful when a slider joint is involved. Though kinematically the two are equivalent 
to a single vector, at times the results obtained by treating them so are not of 
much physical significance. This will be clear from Fig. 23. 

The linkage is represented as a closed polygon of n vectors, with the head 
of one vector at the tail of its following vector. The vectors are termed links and 
the common point of 2 consecutive vectors, a joint. The length of each vector if 
constant, represents Vna link length, and if variable, is a measure of the 
displacement of the slider which may be at one of its ends. There will always be a 
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Fi^ Z.Z. %eior repreterddion ^ linKcLj^es. 


Fi^ 2*5. S^presenidion cf a slider tint 



c. shlual 


(a)l rector 



vector in the direction of reciprocation of a sliding joint. 

The links are numbered starting from 0, which is the input link, through 
n—1. ( n s the number of links in the loop. ) Link n-1 or n-2 is the output link 
depending upon whether or not a fixed link (ie. a vector with fixed length and 
fixed orientation.) is present in the linkage. The fixed link if present, is numbered 
as n - 1. Joints i are the joints at the tail of the i th vector. 0; is the angle made 
by the link i < i.e. the directed vector from joint i to i+1 ) with the positive 
direction of the reference line. (e.g. X axis if linkage is in the X-Y plane) . is 
the angle traversed by the extended link i-1 rotated in anticlockwise direction ( 
positive direction of the axis perpendicular to the plar» of the linkage ) about the 
joint i to get aligned along the link i. For all such cases, the following relation 
holds: 

e,- a e,_i + (2D) 

The angle is constant for sliding and immovable joints and variable for rotating 
joints. It essentially is the the angle between the links i and i-1. 

Number synthesis. 

This involves the selection of a particular combination of revolute and 
prismatic joints from different possible structural configurations. Primary 
consideration is that of the nature of the input and output motion l^^ch decides 
the nature of the first and last joint of the linkage, viz., rotational or 
translational. The nature of the input decides the type of cam-follower mechanism 
to be used, i.e. whether an oscillating or translating follower needs to be 
employed. Oscillating followers generally prove to be better due to the sinplicity 
of ojeration and manufacture. The remaining joints can te selected depending on 


the application. 
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Dimensional synthesis. 


The loop closure equations for the vector representation of th« linkages 
can be written as : 

2 L- = 0. <2.i> 

/ = 0 

where Ly is one of the vectors, and 0, a null vector. 

In terms of the scalar components, 


§ 

/so 

ly COSOy) = 0, 

(22a) 

n-1 

/ = 0 

ly sinCOy) = 0 

(22b) 


where 1; = |L;|, and 6y is the- orientation of ts,- from the reference axis in the plane 
of the linkage. 

The two extremities of the output motion are the two points at which the 
input and output is coordinated. There can be additional points, the total number 
being equal to the number of design parameters. Any number of parameters (Pgfgj) 
which can be either ly, ©y, or ^y can be designed, (generally 2 or 3,) depending upon 
the linkage being used. 


Let Py^ be the output parameter, p^ be the input parameter, and Psec 
the secondary parameter which determines the motion of an intermediate link of 
the linkage. Let P^yt . p^^ . (j = 1,. , . m) be the coordinated input-output 

if 

points .(m is the number of design variables). Then the loop closure equations at 
these points are represented as: 


55 1/ COS(ey) 

/ = 0 





0 , 


(23a> 


and 


[. 




ly Since,) 


^out/ 




]■ 


(2.3b) 
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for 3 B 1, . . . m. 

These are 2m equations in 2m unknowns ( .> and Psec^- ^ which can be solved to 

get the values of these unknowns. The actual equations which are formed fr-om the 
substitution of the above constants (Pq^^ . and Pj^^ .) are dependent on the nature of 
the input, output and the secondarg variable. However, the general method of 
solution remains the same irrespective of these and can be easily taken care of in 
the computational phase. (Section 32 >. 

Note : With a slider joint, it is possible to design the angle also, if needed, and 
can be incorporated in equations (2.3) by means of equation (2T)). 

Inverse analysis of the linkage is done by analysing the linkage formed by 
reversing the sequence of joints in the designed linkage. Analysis of linkages is 
discussed in section 2 £. 

2.42 Caro synthesis. 

Objective of the cam synthesis is to determine the follower parameters and 
the cam profile which will generate a programmed motion serving as the input 
motion to tPe linkage. Design of the follower parameters is dealt with in chapters 
5 and 6. Ciam profile synthesis assLBning that tte follower parameters are already 
designed, is discussed here. 

Systems of coordinates. 

(^m desi^ is done in a temporarily fixed frame of reference, the final 
positicrj of which is decided during the assembly. Here, for convienience, we 
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Fig 2.4: Caw design’ System oj C(yor(lindes. 





Fi^Z.5. Cam syidhesis as vedor sdiillon. 


Finei cam. 



•Vf.: ofthc^otuL reUlire wtceU^ 

Ify jMowcr toting cxniUei 

fttJUtij inside the eatn. 


2-c. Condiiion conUd 



assume that this temporary frame is the global frame of referenoe as far as the 
cam is concerned. 

The follower and the cam are represented in two different frames of 
reference, and the nation of one is mapped onto another. <Fig. 2.4). Cg (same as Tf) 
is the global coordinate system stationed at the cam centre. Lf is the local 
coordinate system at the follower and moving with it. C is the coordinate system 
fixed to the cam and rotating with it. F is the coordinate system fixed in sfjaoe 
and used to excw^ss the follower displacements ( same as ). 


Cam prof ile synthesis. 

Consider cam and follower surfaces in the systems of coordinates described 
above as shown in Fig. 2.4. The mapping of motion from one frame to another is 
done on the basis of the following principle. 

"The position vectors at the point of contact expressed in the global 
coordinate system are equal on the cam and follower surfaces. " 

This also amounts to a single addition of vectors as shown in Fig. 25. 
Along with this, the following condition of contact £23 needs to be satisfied. 

'The relative velocity of the follower with respect to the cam at the point 
of contact should be orthogonal to the common normal to the cam and follower 
surface at the point of contact. " 

This will be clear from Fig. 2.6. 

Using these two principles, the cam profile can be s^rithesised as follows: 
For the point of contact, the globally expressed position vectors are equal. 

?9C = 5gf 


(2.4) 
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where 


and 


?9C 


CR), 


'gc ^ 

are the position vectors of the point of contact expressed in the global franw 
from the oam coordinate system and the follower coordinate system respectively. ( 
represents a vector, and 'C )' represent column vectors. 


CR), 



mge 

= 

[ Me 3 , {R)i^ 

(2.5) 

and 

{R3^ 

S 

[ M/ 3 , {R32f . 

(2£) 

where 





' 

CR>lf 

s 

position vector of the contact pwint in frame L/, 


tM,], 

s 

transformation matrix from frame to C^. 




s 

position vector of the contact point in frame C, 


[M^lt 

S 

transformation matrix from frame C to C^. 


Thus, 

[ Me ] , 

CR>io 

= [ M/ 3 9 (R}2f . 

(2.7) 



{Rile * [ Me 1 [ M/ 3 , {R32f . 

(2.8) 


The above equation gives the position coordinates (x, y, z) of a point on the cam 
profile in terms of the follower displacement and velocity. 


We assume that the global coordinate system is fixed at the cam centre. 
The cam local coordinate system C has its origin at the cam centre and is fixed to 
the cam. r n 


cosOe — sinSc 0 0 

sinSc cos0« 0 0 


t Me ] , 


0 

0 


0 

0 


1 

0 


0 

1 


(2.9) 
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where 6; is the cam rotation angle, i.e. the angle by which C is displaced from Cq 
Thus, 


[ Me ] = 


cos©f sinOf 0 

— sinOf cosdf 0 

0 0 1 

0 0 0 


0 

0 

0 

1 


( 2 . 10 ) 


It may be noted here that on the LJ-I.S. of equation (2.7), [ Me ] » is dependent on 
whereas the R>I.S depends upon the follower displacement and velocity, and 
consequently on the type of the follower being used. [ M# 3j and {R) 2 f 
determined for various types of followers in awiendix B. 


Determination of cam speed. 


The angular speed at which the cam rotates, is determirred by the total cycle 
time of ttie end effector motion. Let T be the cycle time. Then, the angular speed 
is 


2n Nf 

“ = -T- 

where Nf is th« number of cam rotations for one cycle of motion. (Usually, Nf 


( 2 . 11 ) 
= 1 ) 


25 Assembly of motion. 

The following are the objectives of the assembly. 

1. Combination of the motions generated by each CML at the output end. 

2. Coordination of the motions generated by each cam at the input end. 

3. Determination of the assembly dependent parameters. 


Motion combination is achieved through elements known as assemblers, and is 
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essentially a hardware problem. This is discussed in the section 4.1.2 . 


Motion coordination. 

Since the motion in each direction is generated independently, if all the 
motion components at any instant are not appropriately matched, all the Cf’Cs will 
be generating the desired motion individually, but there will be a lack of 
coordination, and the trajectory will appear distorted. This can be controlled by 
maintaining the appropriate phase angle between two cams. (Sect. 4.1.3) 

Determination of assembly dependent parameters. 

Once all the cams in the system are designed, it becomes possible to Icwate 
the cam centre, or the cam axis. All the temporary frames of reference can be then 
transported to appropriate points on the camshaft axis. 

The location of the cam axis depends upon the parameters designed in the 
earlier stages of design, typically, the cam-follower pivot distance, or the offset. 
Generally, the global coordinates of 0^^ (Fig. 2.1) will be known, whereas the 
location of the cam centre Oc^. will be known in terms of 0/^.. This information for 
each cam needs to be coordinated since Vrs global coordinates of all tt% cam 
centres should be such that they lie on a straight line. 

Since only two points in space decide a straight line, this process of 
location of the centres has to be carried out with two cams only. For mechanisms 
with more than two cams, an adjustment in the values of the centre distance etc. 
has to be made in order to maintain a consistency between the mechanisms, l^lith 
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two planar disc cam mechanisms, this problem reduces to a simple geometric 
problem of location of the taints of intersection of a circle and a straight line, 
or two circles, etc. depending upon the type of the follower that is being used. 
(Appendix C.) 

A point to be noted here is that when the origin of tte cam coordinate 
system is transported to some other point by either rotation about or translation 
from the point 0/^., since the position of the linkage remains unchanged, the 
absolute values of the follower displacements may get altered. To take care of 
this, the parameters and 1/ (Fig. 2.7, 2.8) are assigned values at this stage. The 
analytical expressions for the determination of these parameters depend upon the 
location of the cam centre .(Appendix C). 

Zj 6 Analysis. 

Following are the two objectives of the analysis: 

1. To check the accLS'acy of the motion produced by the generated mechanism. 

2. To test the performarwe of the system in terms of 

(i) the peak values of acceleration reached for each member in the 
mechanism. 

(ii) cyclic transmission efficiency. 

The following three step procedure is employed. 

1. Determination at any given instant of the follower kinematics of each Oft.. 

2. Determination of the output kinematics of the linkage at any given instant. 

3. C^jmbination of motion (displacement) in individual directions to determine 

the path followed by the end effector . 
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Determination of the foUoMer kinematics. 


From equation 2.7, Me knoM that 

[ M, ] s {R}j^ = [ Me 3 , (2.12) 

The RJH.S. of this equation is known since [ » f(0f), and ©c depends on the 

lime instant under consideration. AIsud, the cam profile being already designed, 
£R> 2 q is known. The LH.S. is in terms of the follower displacement (s) and velocity 
(s). The above matrix equation is formed of two independent scalar equations with 
two unknowns s, and s and can be solved for the two unknowns. 

Actually, this analysis amounts to solving the loop closure equatins of a 4 
vector loop and can be achieved by treating it as a linkage and using the principle 
for linkage analysis, explained in the next section. 

Linkage analysis. 

The loop closure equations for a linkage are written as 

2 1/ coste,) = 0, (2.i3a) 

; s 0 

and 

X 1,- ElnO;) = D. (213bl 

/so 

With all the dimensions of the lirAcage being designed, for a given value of p^ the 
above two form simultaneous equations in 2 variables, and Pscc^ and can be 

solved using the same techwiique as for the synthesis with the hrelp of equation 2£>. 
Differentiating equations 2.13 with respect to time, we gel lh« loop closure 
equatictfi for the velocity as: 

2 i/ cos(6,) - h sin(e,)-©,- = 0, <2.i4a) 

; 9 0 
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and 

2: 1; sirKe,) + 1,- cos<ei)-e; = 0 (2.14b) 

4 = 0 

Also differentiating equation 2.0, 

e,- . e.-.i + 4,-. (2.15) 

Equations (2.14) are linear in nature, since all the Q/S are determined previously 
form the displacement analysis, and for a specified value of can be solved for 
Potrt Psec • 

Similarly, differentiating equations 2.14 and 2.15 with respect to time, we 
get the loop closure equations for the acceleration as: 

1,- 005(6;) - 2 i; Sin(e;)-e; - 1; COS(6;)e;* - 1; Sin(6;»; = 0, 

(2.16a) 

1; sin(9;) + 2 i; COS(6;)*0; — 1; Sin(6;)6;® — 1; COS(6;)‘6; « 0, 

(2.16b) 

6; - e;.l + . (2.17) 

These are also linear equations, and for a specified value of p^, can be solved 
for Pg^rt and p,*o • 


z 

4 = 0 


and 


z 
<• = 0 


and 
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chapter 3 . COMPUTATIONAL ASPECTS. 


In the previcxjs chapter, a scheme was developed to desisrt the Cl*fl-S in 
different directions as per the motion components of the trajectory. Various 
aspects of this procedure are discussed in the present chapter. These are: 

-Method of solving simultaneous equations. 

-Details of Linkage analysis. 

-Details of Linkage synthesis. 

-Cam profile design details. 

By far, the most common method for solving nonlinear simultaneous equations is the 
Newton-Raphson’s method. A few aspects of this are discussed in appendix A, and 
the relevant issues are refered to in the following sections. Finally, two examples 
illustrating the application of this procedure are presented. 

/ 

3J. Linkage analysis. 

The underlying principle of lirdcage analysis was discussed in section 2.6. 
The Newton-RaiurfTSon’s method is applied to 2 equations in 2 LHTknoMTS, and 
comprises of the following steps. 

1. Identification of the variables. 

2. Selection of an initial point.. 

3. Formation of approximated linear equations in terms of the coefficients of 
the variables and the constants in the equations. 

4. Solution of the linear equations followed by updating the variables, and 
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reformulation till a satisfactory solution of the nonlinear equations is 
reached. 


Identification of variables. 


The loop closure equations for the displacement are represented as: 

ft-i 

5^ cos<©^) = 0, (3.1a) 

; = 0 

and 


ft-i 


; = 0 


1; Sin(e;) 


0 . 


<3.ib) 


Alsoj 


e,- = ©,-_i + 4>j. i = 0, . . . . , n -1 (32) 

The structural parameters of the mechanism are the lengths 1/ of the links and the 
angle between two successive links (if constant). Variables of motion are the 
lengths 1; of the sliding lirdcs and orientations ©y of the rotating links. These are 
dependent of the joint configuration of the linkage. There will be two variables 
which can be evaluated with the two loop closure equations for a specified input 
value. When vectors with varying lengths, and orientations are encountered, it is 
helpful to determine length, rather than the orientation from the loop closure 
equation, and then determine the orientation from equatioi (32). 


An initial point. 

Son^ initial approximate values need to be specified to start tt« procedLff^ 

for the solution of simultaneous equations. For sirwlicity, one of the points at 

which the input-output motion is coordinated during the synthesis, may be stored 

ard be used as the initial point. For analysis over the complete cycle, the 

configuration obtained for every small interval may be stored, and given as an 

initial traint for the determination of the next configuration. This saves the 
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number of iterations to be performed before the solution is reached. 


Fomtaiion of linear equations. 


At any stage of the solution, each term of the coefficient matrix is given 
as: (Appendix A) 


S 



(3.3) 


where e/ is the expression for equation i, 
and Vj is the j th variable of analysis. 

There will be two possible types of values for vj. It will either be the length Ir of 
the link r or the orientation 6r of the the link r. 

If v^- = Ir, 


Sei 

" air ~ 

and C 2 ^- = » 

from equations (3.1). 

Similarly, if v^ = Qr. then from equations 


cosOf, 


sin^rj 


(3.1) and (32), 


(3.4a) 


(3.4b) 




aer 


- i sirrfe,. +P 

p = 0 


(3.5a) 


and Ca^- 


aea 

ddr 


- 2 cos(er+p) 

p = 0 


(35b) 


where z s the number of successive constant ^ps following ©r till r + p < n, the 
number of vectors in the loop, and ©r+F are determined from equation 32. 

Similarly, the constants of the linear equations are given as: 
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G.Sa) 


ft — 1 

bi = — Ip cos(©p) , 

F = 0 

ft— I 

and b 2 a - 2; Ip sirKep) . <3 .6b) 

p c 0 

All the bys and c^-^s are evaluated at cxirrent values of Qo, . , . , . >1»— i), <6 oj - > • » • 
and • I ■ j • .^ft-i). 

Sdlutj43n of equations. 

Tte constants b^- will not be zero unless the solution is reached, and act as 
corrections to the values of the variables V;. 

Thus, 

'^’'next * '^''current ^ 

for the next iteration, where b^- is as evaluated at the end of the solution of the 
linear approximate equations. The condition that provides a check for 

the termination of the procedLa^. 

Inconsistent or divergent solutions reflect that to loop closure with the 
specified values can be struck at. This usually happens when the value of the 
input variable goes out of range of its possible values giving imaginary solutions. 

A few more aspects of the analysis appear in appendix D. 

8»2,Linicage synthesis. 

The principle on which linkage synthesis is based, was explained in section 
2.4.1. To solve the simultaneous equations, the Newton-Raphson's method is employed 
with 2m equations in 2m unknowns, (m s the ntarfser of desigi variables) The 
solution is identified wi%n all the m loop closiv'es are simultaneously satisfied. 
The four step procedure described earlier for analysis holds. 
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Identification of the variables. 


The first fn variables will be the variables to be synthesised. The remaining 
m variables will be the variables Psec representing the tnotion of the 
intermediate link for each configtr'ation j. This also is one of the variables in the 
analysis, and can be found from the method described in the previous section. 

Initial point. 

Initial values are the approximate values close to which the solution is 
expected to lie. In most of the cases, these are formed on the basis of judgement, a 
scaled drawing, or some stw't of approximation. The initial values of the secondary 
variables may also be specified in a similar mamer. Even if arbitrary values are 
assigned, rrast of the time, solutions will be obtained on the basis of the initial 
values of the other variables. However, at times, inconsistencies or divergent 
solutions may be obtained. Usually, with a little trial and error, it is possible to 
reach initial values giving satisfactory results. 

Formation of linear equations. 

The primary variables can either be lengths Ir of the links r, or the angle 
between two lirJ<s r and (r - 1). If the design parameter p^^^j = Ir, the terms in 
the coefficient matrix will be: 


Set 

ai. 


coser; 


(3Ba) 



sin©r; 


Gh) 
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^ " siner; OJd) 

where i - L, m represents the numbers of configur'ations, and <k b 2i — i, 

Jfy-l = ZD, the number of corresponding loop closure equations. 


Similarly, if = i^r, then from equations (3.1) and (32), 


o„^- 




- S sin(e,+p) 

^ tt 0 


(3.3a) 


8e, 


i±i 


2 li'+p cos(e,+p 

p = 0 


(3.9b) 


where z h the number of successive constant ^ps following &/■ till r * p <' n, and 
are determined from equation (32). 


Similarly, the constants of the linear equations are given as: 

ft 

= - 2 cosOp;) , 

P = 1 


(3.10a) 


and — ~ 2 sirk(ap/) • C3.i0b) 

p * i 

All the bfcS and ct/s are evaluated at current values of Qq, . , • , . * 


,%-i)i, and (^ 0 , . , . , . ,^n— i)* for ihe corresponding configurations i. 


Scdution of equations. 

The equations are solved in a manner similar to that described for analysis. 
Inconsistent solutions here reflect that with the specified initial values and / or 
coordinated input output values, to orosistent set of values of the design 
variables can be arrived at. In such cases, either input output crordination, or 
the initial values need to be modified. 


Discussion. 

There are chances of the ooourance of multiple solutions for stfithesis, 
frrare so than analysis, since a larger nun^er of varia bles are involved. However, it 
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usually works out to our advantage that solutions in the vicinity of the initial 
values are obtained. Different initial points may be tried to get different possible 
solutions, and the best one amongst them may be selected depending upon the 
application. Generally, the criteria of space, weight, maximLm allowable link lengths, 
etc. will be the deciding factors. Also, the peak values of the acceleration, 
velroity, and transmission angle will be important. 

3.3 Cam design. 

Cam profile synthesis theoretically comprises of generating infinite points 
that constitute the profile. In practice, a large number of points sufficient enough 
to describe the details of the surface are used. The nimiber depends primarily 
UTKin the the conv>lexity of the profile, and secondly on the size of the cam. For 
most of the cases, 100 to 150 divisions of the cam sixf^ace are sufficient. 
Localised clustering of more points in order to capture the details may be 
necessary in some of the cases. 

Since the profile synthesis equations turn out to be explicit equations of x 
and y in terms of the parameter 6e, <Ref. section 2 . 42 ) these are determined in the 
form of analytical expressions, and quoted in appendix B . 

Analysis of tt% follower motion can be based on a principle similar to 
linkage analysis descrii^d above. 

3J5 ^jplication to a secpnented trajectory of straight lines. <Exaiiple 11 

One of the advantages of using the tase nechanical robots is that non- 
smooth trajectories can be very easily generated. The present example relates to 
a trajectory comprising of a segmented path of straight lines connected to each 
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other. This is shown in Fig. 3.1. The right angled turns present on it, make it 
impossible to be generated with a single cam drive. However, with two CHLs, one 
governing the motion in the X axis, and the oth«r in the V axis, the trajectory gets 
resolved into simple periods of rise, return and dwell for the cam in each CML. 

The resolved components of the trajectory are shown in Fig.s 32, and 3.3. 
The time variation in each section of the motion is assumed to be cycloidal, 
resulting into smooth acceleration ctr^/es. The range of the motion in each 
direction shows that if a simple cam be used to generate these notions, it would 
require an enormous one. This will clarify the use of a linkage. 

The linkage used for motion transmission in the Y direction is an RRF?P 
mechanism. In most of the applications, these will be used on account of their 
simplicity, and the fact that it is iw^erable to have an oscillating follower to the 
cam. The input link of the linkage may be a crar^ or a rocicer. In the X direction, 
an RRPP linkage with a movable support tv^se of assendsler (Sect. 4.12) is used. 
The output displacement simply becomes a sine component of the input lirdc length 
(Fig. 3.4). 

The input to the linkage synthesis is shown in tables 32 and 3.4 for x and y 
directions. The range of the input motion is based on the judgement of the follower 
motion range in that direction. The resulting linkages are shown in Figs 3.4 and 
35. The inverse analysis of these lirricages is done based on the output required 
from the linkages, and is shorn in Figs 3£ and 3.8. This forms Uie input to the 
linkages, and the motion of the follower. 

Both the cams are of the oscillating follower type, ar«J the follower 
parameters have been optimized using the methcai explained in chapters 5 and 6. 
The values of these parameters are shown in tables 3£ and 3.5. Cam profile is 
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based on these values, and the profile generated for both the directions is shown 
in Fig 5 3.7 and 3.9. It is seen that the cam size is not extraordinarily large. 

Analysis of the mechanism in both the directions is done, and the results 
are presented in Figs 3.1D through 3.15. 

3JB Application to a continuous curve. (Example 2) 

The CML type mechanical robots form an ideal choice when the motion in two 
directions exists inherently in a resolved form. It is shown in the present section 
how it becomes possible to generate a series of continuous curves known as the 
Lissajous figures. The principle can be applied to various types of time variations 
of the output displacement. 

Consider the harmonic variation x » a sin(wt), where o is the frequency of 
the wave. If this wave is coupled with a wave in the y direction having a frequency 
which is some multiple of that in the x direction, viz., y = b sinirwt), then these 
constitute the Lissajous figure as shown in Fig. 3.16. In the present example, r is 
taken as 2. 

The variations of the end effectw' motion in the x and y directions are 
shown in Fig.s 3.17 and 3.18 along with the amplitudes and frequencies that were 
used. It can again be seen that the ranges of displacements in both, the x and y 
directions make it infeasible to use a simple cam mechanism. In fact, if just a 
simple cam mechanism be used for the motion in the y direction, one of the 
possible smooth profiles is showi in Fig. 3.19. However, the base circle radius of 
this cam turns out to be excessively large. Fig. 3.20 shows a possible reduction in 
the size through ^stimization, but is still not practicable. A similar phenomenon 
may be observed in the x direction also. The solution to this is the use of 
i 4i 



link'& 9 SSi since relatively larger ditnensions of the links are practicable. 

Tables 3.8 and 3.9 show the iryxit and output motion ranges to be 
coordinated by the lirdcages in the x and y directions. The resulting linkages are 
presented in Figs 321 and 322. A reduction in the lirrfi: lengths may be possible, 
with different combinations of the input motion range and mean value. However, 
generally with large displacements at the output end, large link lengths are 
unavoidable. 

The assembler that is used here is of the x-y projection type (Sect. 4.12) 
which makes it possible to have an RRRP joint combination for both x and y 
directions. Cams for both the directions are of the oscilating roller follower type 
and are presented along with the displacement diagrams in Figs 323 through 326. 

Analysis of motion ftx' the mechanism, and profile analysis of the cam 
surface to find the pressure angle and radius of ctr^ature distribution were 
done, and are presented in Fig. 327 through 3.32. 
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chapter 4. PRACTICAL CONSIDERATIONS 


A conceptual development of the mechanism proposed in the present work 
was shown in the previous chapters. There are other aspects to the practical 
design of these mechanisms. A few of these are described in the present chapter. 

4J. Considerations in the motion assembly. 

In the present section, different types of assemblers, which combine motion 
in different directions, have been dealt with. Also, the effect of the phase 
difference § on the output motion coordination is discussed. 

4A2 Motion assemblers. 

An important requirement of the motion assembler is that it should maintain 
the independency of the motion components along each direction, i.e., the motion 
along one direction should not affect that in any other direction. This is possible 
only if the motion along one direction does not alter the configuration of the 
linkage, or break the cam-follower contact governing the motion along tt® other 
direction. There can be various ways in which this can be achieved, and it remadns 
a task of generating innovative mechanisms capable of doing it. A few possibilities 
of these are discussed here. 

Movable support type assemblers. 

In these, the output end of one of the linkages is rrot grounded. Instead, it 

10592S' 





is mounted on the output end of another linkage. Thus, the assembler serves a dual 
purpose of assembling the motion, and acting as the output link of one of the 
linkages. ( Fig. 4.1 ). This assembler has to be used with a little discretion though, 
since the output end of one linkages moves independently of its input and might 
change the configuration of the linkage. This limits the application of the 
assembler to only those cases where one linkage is of the RRF^ Itfse. Employing an 
RRPP linkage does r»t change the output kinematics due to the change in 
configuration, as can be seen from Fig. 4.1. 

Projection type assemblers. 

Any point on the two dimensional plane can be expressed by the projections 
of its position vector on the two coordinate axes, viz. the abscissa and the 
ordinate. A similar principle can be used with an inverted L shaped assembler as 
shown in Fig. 42. This type of assembler has no restriction on the joints of the 
linkage so far as the output joint of both the linkages is a prismatic joint. 

Splined shaft type assemblers. 

These are used with 2| degrees— of —freedom motion and are explained in 
section 42.1. 

Cable type motion assemblers. 

In these, the end effector is momted on the output «Td of one of the 
linkages, and its motion in the other direction is governed by means of a cable 
connected to the output end of the other linkage. Springs may te used to maintain 
tension in the cable. This assembler can be used in mechanisms generating a 
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planar motion with the polar variables r and 6. Thus is illustrated in Fig. 4.3. The 
independency of the motions is maintained Ur^ough the fact that for small 
rotations of the angle ©, the twist produced in the cable is small enough to keep 
its length unchanged. This type of assembler cannot be used for large angular 
displacements, if accuracy is to be maintained, and in fact is useless when 
continuously increasing output is needed. 

4.13 Motion coordination. 

It was mentioned earlier that the phase angle of each cam governs the 
coordination of the motion between two Ofl_s. This angle is defined as tt« angle of 
cam rotation at time, t = 0. ( denoted by Then, the cam angle of rotation at any 
time instant t can be given as: 

©c = ©«o + wt. 

where u is the rotational speed of the cam. The phase difference between two cams 
p and Q is given by: 

ipf = 

where the notation implies that the cam p leads cam q by an angle 
perfect motion coordination, 

fpf = 0- 

There will be distortions in the the trajectory if this angle deviates from zero. 

It is an interesting observation to see the effect of the variation of this 
angle on the trajectory. It may become useful when the existing trajectory has to 
be slightly changed, or some dwell periods have to te added or reduced. These 
can be easily achieved by manipulating the phase difference. The spatial 
positions of the end effector change drastically with larger values of the phase 
difference. Fig 4.4 shows the variation in the trajectory of example jx^lem i in 















chapter 3 ( Fig. 3.1 ) Miih the phase difference i = iyx ■ It is seen that the shapes 
generated turn out to be very much different from the original shape for which 
the cams were designed, and from each other. 

It may be noted that these different shapes of the trajectory are possible 
from the designed cam profiles without any change in the kinematic characteristics 
( like the cycle time, peak values of the aoceleratiai, velocity, and the extremes of 
the pressure angle, etc. ) of the mechanisms governing the motion in either 
direction. This fact may be made use of to extract some flexibility from the system 
in view of changes in the trajectory. The cbMell characteristics can also be 
similarly manipulated. For example, there were no cb^ell periods in the original 
trajectory, whereas with the triangular trajectory in Fig 4.4, < § = 60, IKD ), at all 
the three vertices of the triangle, there octxff^ a dwell of 2 seconds each. 

Whatever be the phase difference, it needs to be accurately maintained to 
ensure adherence to the desired performarce. 

4.2 Generation of higher degrees-of-freedom motion. 

Most of the discussion in the earlier chapters was related to two degree- 
of -freedom motion systems. In the present section, the possibility of generation of 
motion with higher degrees-of-freedom, based on the same working principle is 
explored. 

A2.± Two and half dimensional motion. 

By 2| D motion, we mean the motion i^Tioh oan be expressed in terms of two 
planar variables, but is 3 dimensional in nature. Essentially, it refers to a motion 
on the surface of a circular cylinder, achievable by the combination of angular 
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motion with linear motion along its axis. Truly speaking, this motion consists of 
only two degrees-of-freedom, but has been treated separately on accoLrtt of its 3 
dimensionality. This type of motion is illustrated in Fig. 4.5. The trajectory is 
produced using two variables - 6 and z, as shown. The mechanism which can 
generate such a trajectory is shown in Fig. 4.6. 

Splined shaft type motion assembler has been used here, and is shown in 
Fig. 4.7. The shaft is hollow, and is free to slide on the axle, and can also undergo 
a rotational motion about it. The splines make the linear motion independent of the 
angular motion. 

Whenever angular motion at the output end is required, a CP/RKP/RXP/R3R 
mechanism may be used for input output ccxjrdination. For small angular motions, an 
oscillating follower may directly be employed. 

It may be noted that the two linkages for the two different directions in 
this case operate in two different planes, and if only disc cams are employed, the 
follower which has to be in the plane perpendicular to the cam profile plane 
cannot be an oscillating follower, since it might cause a loss of contact. A sliding 
follower on the other hand can take motion from the cam in one plane, and can te 
used as the input link of a linkage in a different plane. TTtf'ee dimensional 
(globoidal / cylindrical ) cams can be used in these cases, siroe the plare of 
motion of the follower can be perpendicular to the cam axis. 

Three and more degrees-of-'Preedom. 

Theoretically, the principle on which the Chfl_ ti^e mechanical robots are 
based, can be used to generate motion with any desired nuni^r of degrees-of- 
freedom. In practioe though, it becomes increasingly difficult to confine and 
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coordinate the motion at the output end. The major hindrance is that of the motion 
assemblers. The motion generated in different directions cannot be maintained 
independent of each other with higher degrees of freedom. 

It can be seen that the assemblers bring about the desired independency 
for the applications with two degrees-of-freedom by means of some additional 
constrained motion between two joints which can be considered to be distinct from 
the joints in the four bar mechanisms. ( Fig. 4.1, 42, 43, 4.7. >. This additional 
constrained motion for three degrees-of-freedom, turns out to be a planar motion, 
which is nothing but the projection of the 3-d motion on to a coordinate plare. 
Physical systems exhibiting such a situation can be practicable only if this 
projection is very small. 

If the motion with the third or higher degree-of-freedom be used for the 
purpose of other actuation, like effecting the movement of a gripper, etc. then 
using a third CML might become feasible. 

43 Summary. 

Though a continuous planar motion can be generated using only one CML, 
many a times, the cam profiles so designed are not practicable. Also, non 
differentiable paths are not feasible with these owing to the appearance of cusps 
on the cam profile. 

The major limitation involved in the use of the proposed mechanisms lies in 
the fact that these lade reprogrammabilty, and hence can be used for only those 
trajectories for which they are designed. This Iw flexibility makes them 
inapplicable where the path to be traced by the end effector changes frequently. 
Wherever these changes are limited to a small numba', different cams could be 
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designed with tJTe same linkage for different possible trajectories, and used for 
the corresponding trajectory. Though, changing the cams might become an 
undesirable activity, employing these mechanisms still has its advantages cost 
wise, whenever the number of cams involved is small, and the dtff'ation of the 
operation with one trajectory is large enough. 

The assemblers in some cases may turn out be very cumbersome. The 
dynamics of the assemblers plays an inportant role in determining the performance 
of these mechanisms, and should not be excessively large. This limits the range of 
output motion possible from these mechanisms to as far as the mass of the 
assemblers remains within practicable values. 

It may be noted that these mechanisms are applicable chiefly to those cases 
where a continuous control over the output kinematics is needed throu^Tout . the 
motion cycle. For point to point control, linkages will prove to be better. 
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chapter 5. 

CAM 


OPTIMUM DESIGN OF 
FOLLOWER SYSTEMS. 


5.1 Introduction. 

Profile synthesis of a cam involves the generation of a space curve which 
when traced out by a constrained follower point, produces a desired niotion of the 
follower point. Basic ptrtrase of cam design is fulfilled with the gereration of 
such a curve. However, this does not guarantee an efficient motion transmission 
from the cam to the follower. OU^r operational constraints also need to be 
satisfied for this. These could be classified into kinematic and dynamic 
operational constraints. Two important kinematic constraints are : <i) the kinematic 
efficiency of motion transfer at the higher pair should be as large as possible, 
and (ii) the contact stress at the point of contact should be minimum possible. 

It is seen that both these kinematic operational constraints can be 

measured in terms of critical performance parameters. The pressure angle 

determines the kinematic efficiency, whereas the radius of curvature of the cam 

profile determines the contact stress at the point of contact. Thus the smoothness 

of operation of the cam-follower system can be ensured if the values of these 

performance parameters are withun some critical limits. Host cwnmonly, the maximtan 

value of the pressure angle, and the minimcHn value of the radius of curvature 

depends on the displacement chiaracteristics of the follower . However, i»^ien the 

follower motiwi is maw^ed on the cam profile, structural and othter parameters like 

thie cam-follower pivot-to-pivot distarce, the radius of the roller, thie minimum 
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ihe cam-follower pivot-to-pivol distance, the radius of the roller, the minimum 
displacement of the follower, etc., also play inportant roles in determining these 
values. In the present work, a scheme has been developed in which the values of 
these physical parameters can be selected in a way so as to satisfy the 
operational constraints in an optimum manner, and at the same time, keep the cam 
size as small as possible. 

The current practice in cam design involves checking the values of these 
critical performance parameters with refererce to some roiwgrams, charts, or 
graphs, previously determined for an exhaustive set of values, and has been 
brought out in many publications. Jensen [61, Chen CD, Ganter and Uicker C4] etc 
have explained various ways of providing the necessary checks. Essentially, the 
underlying principle is the selection of a suitable displacement curve and then the 
determination of the range for the best f«rformance, with the help of the charts 
for this displacement curve. These methods have a limitation in that they are 
applicable to a particular curve only, and wt^ the actual ctrve to be traced is 
different, or is a combinatiwi of standard curves, these procedures become 
laborious and inaccurate. Maggiore and Mene^^tti [91 have optimized the cam- 
follower behavioLT by means of a suitable selection of a displacement ctffve 
parameters . Dhande and Kale [21 have developed a method in which a number of 
combinations of design parameters are generated, and finally those are selected 
which yield the minimum size of the cam. 

In the present work, a more formal perspective has been taken and the 
problem is formulated as a classical optimization problem. The cam size is treated 
as the objective function to be minimised with tt® constraints oti the maximum 
pressure angle and minimum radius of cirvature. The ZoutendiSds method of 
feasible directions for nonlinear optimization has been employed (explained later 1. 



The procedure developed is applicable to any displacement curve and has been 
implemented in the form of an interactive computer programme. Only kinematic 
operational constraints have been presently considered. Examples relating to disc 
cam mechanisms are presented. The principle, however, is applicable to other types 
of cam mechanisms as well. 


52 Formulation of the optimization problem of canwfoUoMer systems. 


It is possible to restrict the extremal values of the performance parameters 
within a certain critical limit, but only at the expense of a larger cam size. In 
most of the applications, the designer aims at keeping this to a minimum, since it 
further reduces the weight, balancing problems, space occupied, cost, etc.. Hence, 
the objective ftnction used in the problem is the cam size. By the size of the cam 
is meant the base circle radius. Throughout the work, the terms “base circle 
radius', and 'cam size' are treated as synonymous, since the size of the cam can be 
represented by the base circle radius. 


Let (ai, ag, . . . , aft) be the cam follower system parameters (physical) which 
have to be designed. Since the size of the cam depends on these, let 

r* = f (aj, S.2I ...» Sff) (5.1) 

be the size of the cam. This is our objective function. Let the pressure angle be 
represented as •• 


(ap . . . , a», SjS ) ; 
and the radius of curvature as : 

fi s ^2 (S^lJ . . . , S, S ) , 


(52a) 


(52b) 


where s, k, and s are the follower displacement, and its first and second 
derivatives respectively. 


Now, 
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S = hi<e,) ; i = haOt) ; and s = h#(ef) ; (5 .3) 

where 

&{ = h(t) (5 4) 

is the cam rotation angle. (Usually, ®f = q ; to ■ uniform rotational speed of the 
cam.) 

Thus, eliminatino s, s, and s from eouations 52, 

^ = 3i" (ai, . . . , an, ec ) ; (5.5a) 

P = Sa" <ai, . . . , a„, ©c) . (5.5fa) 

Since we are interested in the extremum values of these furtctions which can be 
obtained as - 


^ dg”i 
def “ dSf 


^ 

ddc d&f 


0 f5.6a) 


d^g»i 

dec2 


< 0 ; 


d2g»>^ 

dec2 


> 0 . 


(5.6b) 


If the value of ©c could be found from these equations in terms of (ai, aa, . . . , an), 
and substituted in equations (5.4), we obtain 

^max = Si'" . . . , an) ; (5.7a) 

Pgjn = 92*” • • • j an) ; (5.7b) 

It is desired that 

i and 1 . (5.8) 

where and are the critical values. 

Let (ajj, . . . , anp and (ai„, . . . , an^) be the lower and upper bounds on the 
parameters (aj, . . . , an) . Then, the problem can be stated as: 

Minimise 

r^ = f (a^i a^* ...» an) (5.9) 
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subject to 


*= 9i'" <«!» . . . , avt) ^ 

^ain ~ S 2 ”* •■■ I a») i 

and 

i *1 i aia 

£ ^2 £ 

a/)| ^ a/} ^ a^^. 

or in a more generalised manner as; 

Minimise 

r^ ® ^ (a^f a^/ . . . > an) 


subject to 


9^' ag) ... I an) i 0> j = j.> . . . j2n+2 

where n = number of variables. 

In words> the problem is stated as follows ; 


(5.10a) 

(5.10b) 


<5. 10c) 


(5.11) 


(5.12) 


"Given the critical values of the pressure angle and the radius of 
ctr^ature, to determine a set of parameters (aju a^, . . . , an) which while 

remaining within given bounds, will keep the pressure angle and radius of 
curvature constraints satisfied, and at the same time yield tt« minimun cam size." 


5.3 Solution of the optimization problem. 

Since the objective furtction as well as the constraints are nwilinear, there 

is no method which can give the optimum solution in a deterministic manrer . What 

instead has to be done is to start from an initial point at which all the 

constraints are satisfied, and go on improving the value of Wie objective function 

in an iterative manner. ThHJS the knowledge of an initial feasible point becomes 
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necessary. It is obvious that the dependency of the objective furotion m the 
various parameters be comprehended before the current values of the parameters 
are altered. In fact, for any family of curves, the direction of the gradient is the 
direction along which the values of the curve constants change at the fastest 
rate. Hence the direction of the negative gradient is always the best direction, 
since we are interested in minimising the furction. 

We may move along this direction as loig as no constraint gets violated, or 
a minima occurs .(Fig 5.1) If a minima occurs, then again the negative gradient 
direction may be determined at that point. If, h»wev&', a constraint is violated, 
this direction becomes infeasible for a further search. At this point then, there 
could be a possibility of finding arother direction along which the objective 
function reduces, may be at a little lesser rate, but without violating any 
constraints. Thus, this direction should be such that 

(i) the objective function reduces along it ; and 

(ii) a finite step in the feasible range can be taken along this direction. 

i.e. it should steer away from the constraint bcHjndary. 

In general there could be several such directions, and we would like to find the 
best one of these. The requirement on this direction - as stated above - is that 
(Fig 52) the Z If should be as small as possible, and Z as large as possible. 
Thus the direction finding problem itself tecomes an optimization problem. 

The step taken along any direotion could be either till a constraint is 

violated, or a minima is reached. It is in^aoticable to directly solve for the point 

of intersection between the line of the direction of search, are! the constraint. 

Hence, in practice, what can be done is to take a step till a minima is reached, and 

then to check if the constraints are violated. If any of the constraints is violated, 

then to recede back arwi modify We step till the constraint bomdary is hit. To 
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■find the step which takes us to a minima, some one dimensional optimization 
technique has to be employed. 

When we see that at a stage either no direction remains along which the 
search can be continued, ex' when the step size reduces to zero, we can assume 
that we have reached the optimal solution. 

5.4 Zoutendijk's method for nonlinear optimization. 

The above discussion forms tt® philosophy of the Zoutendijk’s method of 
feasible directions .(Ref . tiO] for details! In the method of feasible directions, a 
direction is declared as a usable-feasible direction if it satisfies the following 
conditions : 

§T vf (A) ^ 0 

vg/A) 0 ; 

where 

f = objective f metion, 
g^-= constraints, 

§ = direction of search, 

ST = transpose of S, 

A = position vector of the current point, 
n = number of variables, 
m = number of active constraints. 

GsOMtriMlla, this Mans that the angle between the gradient direction of 
the objeotide function and the search direction should be an obtuse angle, since 
only in this direction will the objective funcUon reduce. The second condition 


i = 1, n. (5.13) 

i « i, . . . , n ; 0 = if • - - 

(5.14) 
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specifies a similar condition to ensure that te are always in the feasible range. 
This is illustrated in the Fig. 5J3. 

Going a step further, since we wish to go as close to the negative gradient 
direction as possible, (i.e. as away from the direction of the positive gradient as 
possible) we impose the constraint that 

ST Vf (A) i -a ; (5.15) 

where a is a saclar greater thean zero. 

And similarly since we want to steer away from the constraint to the maximun 
extent, we impose the condition that 

ST vg^-<A) i -Hja ; (5.16) 

where Uj are some scalar multiples of a. Now we see that if we try to irrcrease tte 
value of a, then at its maximum possible value we get the locally best directiqri. 
Then the direction finding problem can be stated as : 

Minimise ( — a ) 
subject to 

ST vf (A) ^ 0 (5.17) 

ST vg/A) - tljOr i, 0 (5.18) 

where (ai, 3 , 2 , ...» an) are constants. 

If this be solved for (sj, Sa, ...» Sn, a), we shall get the optimal direction. An 
important point to note here is that the atove constraints are all linear in nature, 
and so is the objective function. In addition to these, we introduce constraints 
that 

— 1 ^ s,- i 1 » i = 1, . . . , n. (5.19) 

which normalise the search direction. This ensunes that the value of a remains 
finite. Any linear optimization technique can be used to solve this. 

The determination of the step size can be similarly done by any 



dimensional unconstrained minimisation lechruque. After finding the optimal step 
along this direction, it may be iteratively modified if needed to bring the next to 
current point in the feasibility range. 

Again the search is continued till the Kuhn-Tucker conditions are satisfied 
[10] which determine the existereie of an optimal solution of constrained 
optimization problems. In fact it so happens that when these conditions are 
satisfied, the value of a turns out to be zero. (Ref. [10^ Physically this means that 
there no longer remains any direction in which the search can be continued. Hence 
the procedure is terminated. The algorithm of the Zoutendijk's method is presented 
in the form of a block diagram in fig. 5.4 . 

5.5 ^splication of the ZoutendUds procedire to the cam-folloMer optimization 
problem 

The application dependent issues of the Zoutendijc’s method are as follows : 

1. Selection of an initial feasible point. 

2. Objective function 

- evaluation 

- determination of the gradient. 

3. Constraint 

- evaluation 

- determination of the gradient. 

4. Linear optimization technique used for direction search. 

5 dimensional unconstrained minimisation technique for step size 

determination. 

6. Technique used for the modification of the step size. 

All these are discussed in detail in the following sections. 
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Find the gradients of the 

obj. function and the constraints. 


Select (S£-Vf) as the search 


direction. 


Find the optimum direction 5 by 


linear optimization. 



Terminate with (A) 


Apply one dimensional unconstrained minimizatior 

If 

to find the step size (v ). 


as the solution. 


X_ is (>*) 0 ? 


Find the next optimum soln: 
(A);+1 = (A); + V* S 


Is any constraint violated at this point ? 


Find the constraint which is 
violated to the maximum extent. 


Reduce x* till this constraint becomes active 

If 

Assign the modified value to x 


Fig.5.V. 


Algorithm of the Zoutendijc's method. 

















5^.1 Initial feasible point- An initial feasible point is one which has to be 
selected before the Zoutendijk's procedure is invoked. This should be such that all 
the constraints are critically or otherwise satisfied. In our case this point ie the 
value of each parameter of design should (i) be within the range of bounds 
specified for it; and (ii) evaluate to extremal pressLff^ angle and radius of 
curvature values which are also in the desired limits. 

We know that for any given cyclic follower kinematics (displacement^ velocity 
and acceleration) if the cam size is increased, the maximum value of the pressure 
angle decreases and the minimum value of the radius of ctn^ature increases, dil, 
t6)) Thus if initially the size of the cam is made its maximun possible within the 
given bounds on the parameters, all the constraints should be satisfied. If they 
are not, then obviously they cannot be satisfied for any set of values of 
parameters which gives rise to a lesser cam size. To continue with the 
optimization procedure in such a situation, either the bounds on the values of the 
parameters or the critical values of the pressure angle and radius of curvaU-a^ 
have to be modified. 

552 Objective function. 

The base circle radius is used as the objective function. The analytical 

expression used varies with the t^^e of cam follower. The maximum magnitude of 

the radius vector of the point of contact may also be considered to represent the 

cam size, tfrien the foll»aer displacement values are not very large, however, this 

will give results similar to the raie with base circle radius. In all types of cam- 

follower systems, the base circle radius can be determined in terms of the 

follower parameters by means of analytical expressicxTS ret very conel^x in 

nature. Hence once an analytical expression of the ob^ctive function is 

73 



formulated, evaluating it for a specified set of parameters and determining the 
gradient remain simple tasks. 

553 Constraints. 

As far as the constraints putting botrtds on the values of the parameters 
are concerned, determination of the value and the gradient of each constraint is 
fairly easy. Issues become involved when the pressure angle and radius of 
curvatcr'e constraints come into picture. These can be stated as : 

Si “ <«!* •••*«)»>- i 0 <530a) 

arxi 92 = fi* - (ai, . . . , a*) i 0 (52Db) 

To evaluate the value of the constraint, we need to find and for a given 
set of (ai, . . . , an) on the basis of equations 53. Ute see that the evaluation of 
the gradients itself forms an one dimensional (^timization problem, viz., for a 
given set of parameters to determine (D the location (in terms of the cam rotation 
angle) of the extremum values of the pressure angle and radius of curvattre ; and 
(ii) the values of these at the location found in <i>. Also a check becon«s 
necessary to determine whether the obtained value is a maxima or a minima. 

Broadly speaking, there are two ways in which the search can be uridertaken. 
One is the elimination method in which the uncertainity interval in which the 
solution lies is reduced by the elimination of sub-intervals from the range. Other 
is the interpolation method in which the ob^tive function is approximated as some 
simpler polynomial in the expected range of solution and this polynomial function 
is optimized. This procedure is repealed till a satisfactory minimum is reached. 'In 
most of the methods which can te classified into the above two categories, some 
heuristics or approximations are used to quickly arrive at Ihffi solution. An 

elimination method called as the 'exhaustive search method’ however uses no 
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heuristics. It is a so called 'brute fcrce method' in which the objective -Function is 
evaluated at every point in the interval and the minimton selected. This method has 
been implemented in the present algorithm For reasons explained below. 

All other methods are superior to this method only in the view that they 
involve a much lesser number of evcQuations of the objective function (or its 
derivative) thus saving cotnputational time and effort. However, this advantage may 
be outweighed in this particular case by other computaticsTS involved in th« 
search, e.g. the interpolation of the objective function in the interpolation method, 
or the calculations determining the next range of uncertainity to be eliminated in 
the elimination methods. 

Secondly, all these methods though yield a quick solution, .rever guarantee a 
global extremum unless the objective function is a convex function of the 
independent variable. In most of the cam profiles, this will not be so. We are 
interested in the global extremums only, since otiwrwise the results will be of rra 
significance. The method of exhaustive search guarantees a global extrema 
provided that the step size with which the search is undertaken is small enough. 

Since the interval of the search is finite (viz. 360* of cam rotation), we can 
afford to evaluate the objective function at every interval. In the implementation, 
it was seen that this method produces results quickly erwxtgh, and in fact the 
entire optimization procedure gets over in real time. 

Actually, the exhaustive search method could be used as a 'first guess 
method' to locate approximately the stationary point, and then in the vicinity of 
this value the other techniques could have been used. However, it was again found 
that the exhaustive search methal with a reduced step size produces quick and 
accurate results. 
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Thus to determine the extremum, the following procedure is adopted. 

1. Approximate location of the extremum by the exhaustive search method 
using a relatively larger step size. 

2. Precise location using a reduced step size. <Reduced to approximate! iC/i 
of the earlier step sizeJ 

Now let us come to tt« determination of the gradient of these constraints. 
Again, since the constraint itself exists as a one dimensional optimization problem, 
it is almost impossible to analytically determirw its gradient. Hence again a 
numerical technicpje is adopted. Our purpose is to find the vector 


* ‘ ill- 


San 


<52i) 


Each component of this vector, i.e. < ^ ) is obtained by rxjmerical partial 
differentiation, i.e. by varying a,- by a small airount on either sides of the current 
value, - the other parameters maintained constant - and evaluating the gradient at 
those points. This finite difference technique is fouTd to give the gradient, again 

in real time. 


553 Linear optimization for direction search. 


The direction search problem can be stated as follows : — 

Minimise <— a) 
siiDject to 


Si 


93 , 

8ai 


+ 


9aA 


+ lljCt 


i 0 ; 
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sg. 

8ai 


+ 


S3. 

+ Sft — + M.« 

oa. 


i 0 


(522) 


Si 


§f 

8a^ 


+ 


+ 



a ^ D ; (523) 


Si - i i 0, 

s„ - i i 0 ; <5-24) 


-1 - Si i 0, 


-1 - Sft i 0 ; and (525) 

a i 0 ®-26) 

n = number of parameters, m = numter of active ccsTstraints. 

Thte best method for linear optimization presently is the sirw^lex methai and the 
same has been implemented here. However, to ensure nonrregativity of the variables 
used in the simplex method, new variables are introduced, since the present 
variables can take values from -1 to +1. These new variables are (ti,. ... , in) 

which substitute <si, . . . . , Sn) according to the following rule : 

t; = S/ + i ; i = 1, . . . . , n. (527) 

Then the problem becomes: 


Sai 


'C' 


S-' > 

8a,- 
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+ tft 


) (528) 




9^1 


+ 


93. 

+*i.a i 




/ » l.M 


( 


8a, 


ti 



+ tft 


3£ 

da. 


+ 


a S 



af 

8a; 


«^9) 


ti i 2 


tn i 2 ; (5.30) 

li 2 0 , 


tn i 0 ; areS (5.3i) 

a 10. (5.32) 

This problem can now be very easily solved by the simplex algorithm, to give the 
variables (t^, tn, a ) from which either the optimality can te checked ( a = 0 ) 

or the direction of search obtained (s, = t, — i ). In the implementation, all the 
values 'tlj are taken to be unity to give equal importance to all the constraints. 


5.5.4 One dimensional unconstrained minimisatian for step size determination. 

Let 

r^ ss f (aj,, aj, » an) (5.33) 

be the objective function in 'n' variables (ai, ajj > *«) . 

If S s (si, . . . , s.)T is the direction of search and X is the step size, then the 

next point can be written as : 

ss A, + X S = C (aj+Xs^), (a/j+Xs^) 1^. (5.34) 

and the objective function as ; 


r* = f ( (aj+XSi), 


(an+XSft) ) 


(5.35) 
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where (ai, a 2 . 


t aft) and (sj, . . . , Sn) are constants. 


This is a single variable function in X which needs to be optimised without ajTy 
constraints. 

As previously stated, one dimensional methods can be of the interpolation or 
elimination type. The interpolation methods are based more on mathematical 
approximations of the functions than numerical hetr'istics, and they converge after 
a finite number of steps to the real solution (Ref CiOD. They are also more 
efficient than the elimination methods. Hence it will be better to use one of the 
interpolation methods. Since in our case, there are possibilities of finding one 
optimal solution close to the previous one, a method which gives quick results in 
close vicinities of the current point is necessary. The direct root method is one 
such method. It also works out well for long distance searches though with a little 
lower efficiency. In this method, the derivative of the objective function is 
approximated as a straight line and the point at v^dch it becomes zero is found. 
Then this becomes the new optimal point. This procedLff^ is repeated till a 
satisfactory value is obtained. This is nothing but determination of the root of 
the equation 

Stj 3f( (aj,+Xsj,)) > (an+XSft) ) ^ 

^ = ax " ' 

by the Regula — Falsi technique. 

5J5J5 Step size modification. 

On taking an optimal step along the optimal direction, if any of the 
constraints get violated, then the step size is reduced till the constraint 
boundary is reached. The maximun violated constraints is first considered, and tl^ 



process repeated till all the constraints are finally satisfied. The procedure uses 
a simple linear interpolation between the feasible and the infeasible point; every 
time evaluating the gradient to check whether the value becom es zero. 

If the rth constraint is violated at X ■ X*, 

OV = flf I X = i 0 satisfied at X = XiC = 0, initially) 

(537) 

g"r * ar I ^ = gf(A^ + X*S) i 0 violated at X = X* 

(5.^) 

Then, assuming a linear variation of gr wr.t. X, 

g,(X) = p^ + p^X (5.39) 

where 

gv — g'V 

Pi = gv, and P 2 = — (3.40) 

X* 

The approximate value of X which makes g^CX) = 0 is the desired value of the 
modified step. 

Then, Pi+ PjX = 0 

X= X* (5.41) 

Pj — s'V 

Using this value, a new point can be fotflid which may be in the feasible or 
infeasible region, A similar interpolation as above can then be contimed till the 
next point lies on the constraint boundary, i.e, 

gr (X) 1^ . s: 0 (5.42) 


5 £ Lifnitations of the 2toulendi3k's method . 


Though the ZoutendiSc's methcd appears to work out efficiently, there are 
certain drawbacks. These may be listed as follows : 

1. It involves a lot of computaticn in finding the optimal direction owing to 
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the linear optimization. 

2. It requires the evaluation of the gradient of the constraints and the 
objective function increasing the computation further. 

3. Similar to other nonlinear optimization techniques) a global optimum is not 
guaranteed by the Zoutendifl<’s prooedLff^ and it needs to be repeated with 
different initial points. 

In spite of this, the direct approach used in this method which relies on the 
dependence of the search direction on the nature of the objective fuTction and the 
constraints, the method works out efficiently. 

In the following chapter we shall discuss hwsw it has been applied to typical 
cam-follower mechanisms. 
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chapter €. application to typical 
CAM-FOLLOWERS. 


6JL Rreliminaries. 

In the application of the method develoc^d in the previous chapter to any 
specific cam-follower mechanismj the following issues are important. 

1. Identification of the design variables. 

2. Formulation of the objective ftsTction. 

— (a) Evaluation of the function value. 

— (b) Evaluation of the function gradient. 

— (c) Evaluation of the function derivative for single dimensional 
optimization. 

3. Formulation of the constraints. 

— (d) Evaluation of the constraint value. 

— (e) Evaluation of the constraint gradient. 

The objective function as previously mentioned will always be the radius of the 
base circle of the cam, and the design variables will be those which decide this 
base circle radius. Analytical expressions for all the ttr^e (a), (b), and (c) listed 
above will be used; which will be typical of the cam under consideration. To have a 
generality in the evaluation of the pressure angle and radius of curvature 
constraints, the 'transmission angle' instead of the presstre angle has been used. 
This is nothing but the complement of the pressus^ angle and is ttefined as: 

T BB (^ — 0 ) and (6.1) 




(62) 
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where it is the pressure angle. 


Then the value of tlie pressure angle constraint will be 

^ 0 ; 


(6.3) 


where 


T* = (^ - 0*) 

is the critical value of the transmission angle. 
Similarly, the radius of curvature constraint will be 


fi - A 


Blit 


0. 


(6.4) 


(6.5) 


It may be noted that the critical value of the radius of curvature here refers to 
the positive value of the radius of Cfj^ature. For the negative value of the sans, 
the condition is that the absolute value should te greater than the radius of the 
roller ( for roller followers ), or should never be reached. (Flat faced follcft'iers). 
This may be incorporated as another constraint, or checked with the final solution 
obtained. 


It is always necessary to delimit the range of permissible values that the 
variables would take during any stage of the optimization, for at times the 
variables will take values which will have no physical significance, or impossible 
to arrive at. Hence for each parameter, the upper and the lop«r bounds on the 
values need to be imposed. These act as additional constraints during the 
procedure of optimization; and are evaluated as follows: 



i. 

0 , 

and 

(6£) 

a,- - a,-, 


0. 


(6.7) 


A fact to be noted is that all constraints evaluate to a value lesser than zero 

when satisfied, and to a value greater than zero when violated. The gradients of 

the constraints forming the bourds on the values will obviously be 

Vg; = ( 0, 0, . . , z, . . , 0 (6.8) 
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where 2 will have a value ±1 depending upwi whether the constraint under 
consideration is the upper or lower bound on the value of the variable. 

The gradients of the pressure angle and the radius of curvature will be 
evaluated by the finite difference method as explained below. 

Let C Aai, . . . , Aa^ JT be the steps taken on either side of the CLff'rent 
values C ai, . . . j an respectively, for this finite differentiatiwi. Then for 
a constraint g^- the gradient will be given as: 


g/ai,..,(ay+aa;),..At)-gj-(ai,..,(a,-Aa;),..,a») 


(6.9) 


Each term of this vector demands two evaluations of the constraint, and in all for 
one evaluation of the gradient, 2n evaluations of the constraint will be necessary. 
This means that the pressure angle or the radius of oLrvature needs to be 
evaluated 2n(N + Np) times where N = the number of divisions on the cam surface, 
and Np «= the number of evaluations done for the precise location of the extremum. 
This figure though alarming, in practice turns out to be negligible as will be seen 
from the time taken for one iteration of the complete procedtr^. 


The following sections illustrate this idea in relation to typical cam 
follower mechanisms. In the presentation of the results, the following format fw' 
the input and output has been followed. The input to any problem will be in terms 
of the following: 

1. An initial feasible point. 

2. Upper and lower bounds on the values of each variable. 



3. The critical values of the pressure angle and radius of curvature. 

4. The step size necessary for each parameter to evaluate the gradient of 
the constraints. 

Output of the procedtr^ will be in the form presented in tables 6.4 and 6.6. 
The first column shows the values of the parameters with which the iteration was 
started. This will always be a feasible point, and the values of the constraints at 
this point are shown in the next column. The third column shmws the direction of 
search from this point, and the fourth one, the optimal step that should te taken in 
this direction to find an one dimensional minimtm of the ob^tive fuxjtion. 
However, at the resulting point (showi in the next colurm) the constraints may be 
violated, (Values of the constraints at this point are stwDwn in the next coltann. ) 
and therefore, the step has to be reduced to a value shsjwn in the colLonn 'modified 
step’. The resulting point is shown in the first colutm of the next row, and the 
values of the constraint in the third row. Finally, the last colLsnn shows time 
in seconds at the end of each iteration. 

We shall now consider the various cam-follower mechanisms. 

62 Translating roller follower mechanism. 

In this section the principle of optimisation as applied to a translating 
roller follower mechanism has been explained. (Fig 6.1). 

62.1 Selection of design variables. 

The bese circle radius of a cam is delarmined b« the minimum fcUower 
disolacemenl, the roller radius and the offset. Hence, these are seUcted as the 
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hU-lfarandcrs ^ ihe raiprocaiin^ 
roller follower. 


■<r 



fig 62 Paramekn 

roller /o! lower. 



design variables. In actual practice one of these may be a constant depending upon 
the application of the cam. In such cases these constant values may be assigned to 
the corresponding parameters, and the cam optimized for the remaining ones. 
Typically, the radius of the roller is most of the times decided by constraints 
other than the pressure angle and the radius of curvature, for e.g. - the force 

transmitted, the minimum pin radius, etc. Here all of these have been considered 

for the purpose of generality. 

622 Formulation of the abjective function. 

At the minimum follower displacement, the roller follower will be tangential 
to the base circle of the cam. The prime circle radius will then be 

rp = 4 ®* + ^ 

and the base circle radius 

r^ a 4a* + Sq^ — rr <6.11) 

The gradient can be obtained by differentiating this wr.t. a, Sq and r^. If X is 
the step size along the direction <do, dj, d 2 ) then, 

r^ = 4^+^ + <So+Xd qJ* — fr'f+Xd a) (6.12) 

When differentiated wr.t. X, we get r/(X). 

All the resulting expressions are shown in Table 6.1. The objective function 
reflects that the cam size increases with increasing 'so' and 'a' and decreases with 
increasing 'rr'. However, with a reduced roller radius, the permissible value of tire 
minimum radius of curvature may have to be reduced, ( since (r should te 

always less than rr ) or r* may have to be decreased. This may result in a 
reduced cam size which may violate the presstre angle arrf / or radius of 
curvature constraints. 

The radius of the prime circle could also have been used as tte objective 
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expressions used for TRAfGLATTNG ROLLER FOLLOWER. 


Design variables: 


minimum follower displacement 

So 

follower offset 

a 

roller radius 

Pr. 


Objective -function: 

value: 



gradient: 


So 

“ "I 

12 2 

» 2 ^ 2 

1 

+ 

0 » 

1 

^So + a . J 


single variable function for l_din)ensional optinuzation 


r^CX) = ^ (S(3+XdjjJ + ta+Xdj)^ ~ 

, (So+Xd(j)do + (a+Xdi>di 

r^'<X) = I 2 ' 2 ■ 

(sQ+Xdo) + 


(r'r+^ci2) • 

dg. 



■function, esp@cially in the cases where the roller radius is decided beforehand, or 
when the lower and the upper bounds on the roller radits are very narr'ow. 

G2.3 Pressure angle and radius of curvatu^ constraints. 

The cam follower geometry has been analysed (Appendix B) and the 
expressions for these have been determined analytically It is seen that the 
pressure angle is independent of the roller radius. This is due to the fact that 
the actual point of contact of the cam and the follower i.e. the IQr (Fig 6.1) is 
governed by the offset and the minimLsn follower displacement only. The expression 
for the radius of curvature reflects the complexity involved, and it is seen that 
the nonlinearity of the various terms make it almost impossible to ccwnprehend the 
precise effect of the variation of any ore of the parameters on the minimum 
radius of curvature. 

62.4 Solution of an example problem. 

The input given to the example prchlem demonstrating the optimization of a 
translating roller follower is shown in table 62. Practical values have been 
selected. Gererally, the pressure angle values for such a follower shtxild not 
exceed 30* to 35*. The minimum radius of curvature should be such that the minimum 
prime cur^e radius turrs out to be mor^ than the roller radius. The step sizes 
inputted to the problem are based on judgement. Cfenerally, it turns out that these 
produce satisfactwy results if selected to be about i to 2 % of average 
possible value of the variable in that range. The range of the permissible valiffis 
again depends upon the application. 

The output / results of this are p-'esented in table 6.4. The cam size has 
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step sizes: 


St^ sizes: 


l.Hin, follouer di^lacefient (s 6) 

e.i 

l.Pivot-to-centre dist^ (a) 

0.1 

2. Offset (a) 

0.1 

2.101101^ (l_f} 

6.1 

3. Roller radius (r_r) 

ai 

3. Rolls' radius (r.r) 

6.1 



"lliin follouer di^la«nent (beta) 

6.2 

fioijids on values: 


[kMfe on values: 



1 u 


1 u 

l.Mifi follouer displacepient (s 0) 

3.8 20.0 

l.Pivot-to-centre distoce (a) 

5.6 25.0 

2.0ffset (a) 

0.0 18.0 

2.Follo«er K'ft (l_f) 

5.6 25.6 

3. Roller radius (r_r) 

8.5 <8 

. 3.Roller radius (r_r) 

1.6 4.6 



4. Min fol lower displacement (beta) 

15.6 40.0 

5.Hin. radius of arvatire 

1.0 

5. Mia radiiK of arvature 

1.6 

8.tex. pressir^ ugle 

35.0 

S.Hax. pr^siFeai^le 

37.5 

Initial values. 

1 

Initial values. 


l.Hin. follouer displacement (s 0) 

16.0 

l.Pivot-to-csitre distsice (a) 

16.8 

2. Offset (a) 

5.0 

2. Foil ewer arm (l_f) 

5.6 

3. Roller radius (r_r) 

1.0 

3.Roller radius (r_r) 

1.2 



4.fiiniBiffi follower displacement (teta) 'KI.O 

Note: fill angles inote^. 


Note: fill sigles in degr^. 


Table G.eirpt data for translating follouer. 

T^le 6 . 5 Input data for osclllatii^ follouer. 





87.2922 (-15.2fS^ -22.6W2 82.72^) (34.3890 -3.6212) 

1.3538 ( 18.4387 3.1847 2.9315) (26.9874 5.8517) 

18.9900 ( -1.7425 -8.5328 15.6594) (89.9733 -8.2294) 

4.5324 ( 4.8764 -1.6598 1.1788) (49.1757 2.1958) 

2.4647 ( 8.4939 2.5963 4.7537) (28.7966 1.7338) 

3.4684 , ( 4.4751 -1.6178 1.9233) (59.8392 1.9865) 

2.5816 ( 8.2362 2.8959 5.6644) (39.6424 8.6219) 

2.8365 ( 4.2285 -1.6858 2.3591) (52.8905 8.3563) 

9.6683 ( 5.6822 1.1933 3.9751) (33.8575 1.1867) 

9.8312 (-1.1211 -8.2354 19.9268) (^8975 -1.5931) 

9.5877 ( 5.9331 1.7634 2.7722) (35.3443 1.8381) 

1.5665 ( 6.6398 8.1686 4.3181) (34.9211 8.4782) 

I 

1 


1.4595 2 
1.3538 4 
1.5111 6 
2.7419 9 
1.7199 12 
1.5284 14 

8.9172 16 
8.7964 29 
8.6683 23 
8.8891 26 
8.5535 39 

9.1512 33 

1 

I 


9.5376 ( 4.5364 2.2721 2.6428) (40.5239 ' 1.8729) 8.9906 56 




Ito, Current point Constraints Size Search Direction 
s.zero ofPset r_radius pr.ang. rad.c(jrv. s.zero offset r.r 


9. (19.0098 5.0900 1.0990) (3^ 3736 6.7997) 10.18 (-0.6325 -0.3162 0. 

1. ( 9.0769 4.5385 2.9329) (34.9999 4.9844) 8.12 ( l.OCOO -1.0000 0. 

2. (10.4397 3.1847 2.9315) (26.9674 5.0517) 7.97 (-0.6763 -0.2065 0. 

3. ( 9.4088 2.8726 4.6099) (27.8103 3.1812) 5.84 (-1.0000 -1.0000 -0. 

4. ( 6.6678 0.1317 2.2891) (35.0000 2.5282) 4.38 ( 0.7409 1.0000 1. 

5. ( 7.9355 1.8426 4.8000) (26.6256 2.12G6) 4.15 (-1.0000 -1.0000 -0. 

6. ( 6.4871 9.3142 3.9828) {^.m 1.5243) 3.33 ( 0.7085 l.OCOO 1. 

7. ( 7.9569 1.2315 4.8090) (28.7703 1.3703) 3.16 (-1.0000 -1.0000 -0. 

8. ( 6.3595 0.5250 3.5913) (34.1888 1.0800) 2.78 (-1.0000 1.0000 -0. 

9. ( 5.6822 1.1933 3.9751) (33.0575 1.1067) 2.73 (-0.6920 -0,1453 0. 

19. ( 5.6268 1.1817 3.1317) (33.3138 1.0000) 2.62 (-1.0000 1.0000 -0. 

11. ( 5.9734 1.7S1 2.7932) (^.0000 1.0092) 2.57 ( 1.0000 -1.0000 0. 

l 

' t 

I 

18. ( 5.8749 1.7345 2.8826) (34.3939 1.0000) 2.48 (-1.0000 1.0000 -0. 

19. ( 5.6734 1.7351 2.8823) (35.06®0 l.lKfll) 2.48 

Opiimitation iht translating SolU>mr. 


reduced from 10.18 lo 2.48 in 19 iterations. Many a times this minimum value itself 
might turn out to be so lal^l that it is impracticable. In such cases^ when the size is 
determined by other constraints like the minimiffn shaft diameter etc., the minimum 
permissible value of the size as regards those constraints may be safely asstaned. 

It is seen that the values of the pressure angle and the radius of 
curvature are critically satisfied. In some cases the constraints may be satisfied 
■ critically or otherwise - even at the lowest possible cam size. Such a condition 
should always be checked before the optimization procedcs^ is invoked, since in 
these case it will not be necessary to optimize at all. optimization can, 
however be tried with the constraints tightened w' the boards on the values 
changed. 

It is seen that the direction of search in n»st of the cases is found by the 

linear optimization technique. (Reflected from the fact that This shows 

that the constraints were critically satisfied at that point. It will be seen from 

the step size in the preceeding iteration, that the optimum step size has been 

modified to take into account the constraint value. Wherever no constraint was 

violated, the optimum and modified step sizes are equal. The table shows that trust 

of the time the constraints were violated. The modified step shows them to be 

* 

satisfied. 

A look at the last column shows that the total time taken for the 19 
iterations was 56 seconds. The results were taken on Apollo-Nexus 3CH30 terrrunals, 
which are considered to be moderately fast machines. Thus we see that the large 
number of evaluations of pressure angle and radius of ou^ature, on the i*dTole, 
becone negligible. 

Fig 6.5 shows the variation of the pressure angle for one conplete rotatiwi 
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of the cam designed on these results. The radius of curvature variation areil the 
cam profile are shown in the figures 6.6 and 6.4 respectively. 

6.3 Oscillating roller follower mechanism. 

In this section the optimization procedure as applied to an oscillating 
roller follower mechanism has been explained. Ref Fig. 62. 

63A Selection of the design variables. 

The base circle radius fo the cam here is determined by the cam-follcM*4er 
pivot to pivot distance ^a', the follower arm length 1^', the roller radius "r/ and 
the minimum follower displacement 'iS'. Again, one or more parameters may have 
preassigned values. This can be taken care of in a manner mentioned earlier. 

632 Formulation of the abjective function. 

At the minimum follower displacement, the roller is tangential to the base 
circle. From aAOB we get 

rj** - a® + — 2 a 1/ cosjS (6.13) 

r^ = r^f — rr = -ja^ + 1^ — 2 a 1/ cos$ — rr. (6.14) 
This forms the objective function for this particular tyse of follower. The 
gradient can be obtained by differentiating this wr.t. a, If, 0, and rr. Resulting 
expressions and the expressions for single variable optimization are shtMMn in 
Table 6.5. The cam size is seen to increase with a, If, and 0 (when 0 i ~>. However 
it reduces with the follower radius. 

The prime circle also could have been used as the objective function if the 
roller radius was preassigned, or its range was narrow. If the value of 0 ccxild be 
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T«$le 6S, 

expressions used for oscillating roller follower. 


Design variables: 

pivot-to-cam-centre distance 
follower arm length 
roller radius 

minimum follower displacement 

Objective function: 

value: 

r^ = ^1^ + a— 21/a cosg — rr. 



a — 1 /cos(3 



l/^+ a^ — 21/3 cosiS 

gradient: 


1 /— a cosi3 


N 

2 2 

1/ '+ a *— 21^a cosjS 



- 1 



a 1 /SinjS 



1 2 2 

[1/ + a — 21/a CO3J0 





Single variable function for l_dimensional optimization: 

r^OJ = ^| Q /+Xdo)^ + (a+Xdi)^ — 2Q z+Xd©) (a+Xdi)cos(j3 + Xdg) — (rH-Xd 2 ) 

_ /+Xd(j)do + (a+Xd^ldj— p ^ 

\ a ) — ■ — ■ ■■ ■ ■ 

^ Q ^+XqIo)^ + (a+Xdj^) 



where p = 1 /dg sin(j3 + Xd^) — (d©! / + s. + Xd 3 ). 



fixed up from other considerations, then an expression similar to that of the 
reciprocating follower would have been obtained. 

6.33 Pressure angle and radius of curvature constraints. 

The pressure angle unlike in the previous case is not independent of the 
roller radius. This is due to the fact that the direction of velocity at the point of 
contact on the follower depends on the roller radius of the follower. The 
expressions again are nonlinear in nature, and do not reflect any direct 
conclusions regarding the deperwlenoy of the pressure angle and radius of 
curvature on the design variables. 

63.4 Solution of an example problem. 

Input to the problem is shown in Table 6.3 • Normally, the permissible 
values of the pressure angle for an oscillating follower are more than that for a 
translating follower. Values of j3 are accepted in degrees, but during optimization 
are calculated in radians. The step size for the gradient evaluation for each 
parameter is consistent to the mean value variation. 

The output / results are presented in Table 6.6. Cam size has retfciced from 
5.76 to 2.48 in 54 iterations. This low reduction in size is chiefly cfeje to the 
presence of the parameter $ in the objective furcticwi and as one of the desi^ 
variables. The difference of scales in the values of $ and the other parameters 
restrict the step size to a very low value since a larger step violates the 
constraints. The objective ftffwtion also gets distorted cfcje to these disparities in. 
scale. This could be taken care of by normalising ard rescaling the variables. If 
dimensionless parameters of some sort be used, and all the variables normaUsed, 
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No. Current point Constraints Size Search Direction 

a l_f r_r beta pr.ang. rad.curv. a l_f r_r beta 

e. (18.0889 5.8900 1.2888 8.6981) (37.3375 4.3999) 5.76 (-1.8080 1.0808 -im -1.808 

.1. ( 9.9963 5.8837 1.1963 8.6945) (37.5880 4.3859) 5.74 (-1.8880 1.0898 1.0098 8.334 

2. ( 9.9854 5.0146 1.2873 8.6981) (37.2692 4.3782) 5.73 (-1.0880 -1.0088 -1.9880 -1.800 

3. ( 9.9813 5.0186 1.2033 8.6941) (37.5888 4.3609) 5.71 (-1.8800 1.9008 1.9980 0.333 

4. ( 9.9692 5.0227 1.2154 8.6981) (37.2488 4.3428) 5.71 (-1.8880 -1.0800 -1.9888 -1.888 

5. ( 9.9648 5.0183 1.2189 8.6937) (37.5808 4.3315) 5.69 (-1.8880 1.0900 1.0808 0.333 

6. ( 9.9514 5.0316 1.2243 8.6981) (37.2258 4.3126) 5.68 (-1.8880 -1.0808 -1.0080 -1.880 

7. ( 9.9466 5.8268 1.2194 8.6933) (37.5890 4.3882) 5.66 (-1.8809 1.8880 1.8080 9.332 

8. ( 9.9320 5.8413 1.2340 8.6981) (37.2801 4.2795) 5.65 (-1.8880 -1.9900 -1.8000 -1.800 

36. (9.2^ 5.4874 1.^1 8.6981) (35.9391 3.1051) 4.57 (-1.8900 -l.WO -1.0000 -1.800 

37. ( 9.2074 5.3794 1.5721 8.6701) (37.5800 3.8209) 4.43 (-1.8008 1.8888 1.0080 0.286 

[53.( 8.0582 5.9829 2.1755 8.6892) (37.5899 l.eCfiS) 2.48 (-1.8880 1.0080 im 0.28S 

54. ( 8.0577 5.9834 2.1761 9.6893) (37.4915 1.8080) 2.48 

0/ iht~ osciiUHnf~/oiio'^er' 


a 


Opt.step Next point Constraints Hodstop lifle 
3 IJ r_r beta pr.ang. radcurv. 
e.7627 (3.2373 5.7827 6.^73 -6.66^8) (89.9877 8.8876) 8.8837 19 

8.5722 (9.^42 5.5758 1.7685 9.8856) (28.7586 3.7979) 8.8118 26 

0.5796 (9.^ 8.6277 8.1185) (78.8991 2.2639) 8.06« 52 

8.5717 (9.4)96 5.5823 1.7758 8.88^) (28.7868 3.7738) 8.8121 57 

0.5803 (9.3889 d-M2'l 8.6351 8.1173) (78.9«9 2.2265) 8 . 89 ^ 76 

8.5676 (9.3972 5.5859 1.7786 8.8827) (28.7791 mi ) 8.9133 81 

0.5809 (9.3706 dW 0.6^4 8.1173) (78.97^ 2.1878) 99 

8.5655 (9.3811 5.5923 1.7858 9.8812) (28.7981 3.7215) 8.81« M 

8.5815 (9.3585 4^99 6.6525 9,1166) (79.8999 2.1453) 9.8853 124 


9.6059 (8.6295 4.8915 8.9941 8.8922) (88.3177 1.2154) 8.8288 438 

0.5671 (8.6403 5.9465 2.1392 8.8328) (38.7368 2.4682) 8.8978 435 


0.5377 ( 7.5205 6.5296 2.7132 8 . 7228 ) ( 35.5961 9 . 5032 ) 0.8885 - 528 

m mmww 031 O 










i- Specification of initial values/ bounds and step sizes. 

2. Effects of different displacement curves. 

3. Methods of improving the algorithm. 

6J5.i Specification of the data. 

(D Initial values of the variables. 

As already mentioned, it is necessary to give an initial -feasible solution to 
start the process of optimization. In most of the cases, this will be selected as 
the one yielding the maximum possible cam size. A little trial and error may be 
necessary before these values are finalised. An important point is that the 
solution obtained is always dependent on the initial point specified. However, this 
dependence might not be always predictable. The only alternative is to use 
different initial points and obtain different results if possiblej arel use the best 
one. It is to be noted that the solution obtained will not necessarily be a global 
optimum. 

(ii) Bounds on variables. 

Practical constraints do generally put some bounds on the values of the 
variables. In cases where such tounds are not precisely defined, arbitrary bound 
values may be specified, but the procedure does not distinguish between loose 
bouTKls and tight bounds. i.e. no distinction between flexible and rigid constraints 
is made. All constraints are treated as rigid, and are not allowed to be violated. 
Practical situations may be different, and the fcKxnds may be loosened or tightened 
accordingly. 
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Generally; lower bounds on the variable values will be rigid owing to 
physical limitations. Upper bounds will mostly be placed in view of size reduction, 
and will have to be varied if the constraints ( pressure angle and radius of 
curvature ) are not satisfied. 

(iu) Number of divisions for constraint evaluation. 

The accuracy to which the constraint gets evaluated depends on the number 
of divisions of the cam. Generally, 80 - lO) divisitms of the 360 * at the cam centre 
prove to be a good choice. This depjends on the application, and larger cams may 
demand larger number of divisions. 

(iv) Step size for gradient evaluation. 

In evaluating the gradient by finite differentiation method, the step sizes 
< defined as the vector ( Aa^, . . . , Aan 3 ^ ) piay invDortant role. If the step 

size is too small, the value of that gradient component will turn out to be zero. If 
too large, it will yield incorrect directions. A little judgement in this view is 
necessary tefore these values can be properly decided. Generally a step size of i 
- 2 of the mean value of the variable in the specified range wotH<s well. 

63.2 Consideration of different displacement curves. 

The values of the cons traints obtained at any stage will be to some extent 
dependent on the type of the displacement curve ciSed. HoiAjever, these are not of 
much significance if the second derivative, viz. the the acceleration is continuous. 
Otherwise, the radius of curvature values become very small yielding cusps, and 
sharp corners on the cam profile. Displacement curves having smooth acceleration 
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variation yield the best cam surface characteristics and therefore, the radius of 
curvature. 

655 Methods of improving the algorithm. 

(i) Using heuristics for locationg the extremum values of the pressure angle 
and radius of curvature: This is done in view of reducing the computation during 
the cyclic evaluation of the constraints. One possible heuristic will be to search 
for the extremum in the vicinity of the location at which it occured in the 
previous search, and use some quick search tectnique to locate it precisely. 

(ii) Using other techniques for step modification: Since the point at which 
the line of' the search direction intersects the constraint boundry is found by 
iterative interpolation, it takes much computaticnal time. Using an accelerated step 
technique, or fitting a quadratic might yield a better solution in some cases. 

^ Checking the order of constraint evaluation: If the constraint which 

would be most commonly violated is known, then this constraint should always be 
evaluated first during the modification of the step size. Any other constraint 
should be checked for violation only if this constraint is satisfied. 
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chapter 7. IMPLEMENTATION 


The theory regarding the design of the cam mcM±jlated linkage type 
mechanical robotsi followed by the optimization procedure to minimize the cam size 
was explained in the previous chapters. A software incorirorating these two has 
been developed and implemented on the ApoUo-Nexus 3CBS3 ( b & w ) workstations. 
Details regarding the implecnentation are described in this chapter. 

The system has been developed in C, and primitive graphic resoLff'ces 
provided on the Apollo workstations have been used fw' all the graphics involved. 
The system developed is interactive in nature, and provides user control at all 
stages of design. 

7J. Scope and assunptions. 

The prcM^Lff^ of optimization is included in the system as a part of the 
design of the CMLs. Hot«ver, the optimization procedij^ can also be invcrfced 
independently. 

It is assigned that tt® specified trajectories will always be two variable, two 
dimensional trajectories, in the x-y plane. Though mechanisms for more than 2 
dimensional trajectories can be designed, they will not be displayed in three 
dimensions. Instead, they will appear as x-y plots in 2 dirrensions. 

The cams that can be used are of the disc type, with oscillating and 

reciprocating, and roller and flat faced follower combinations. Four bar linkages 

( prismatic or revolute joints ) with the number of vectors upto eight forming the 
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loop can be designed. ( More than eight vectors will almost never be needed. > 
Optimization of the cam-follower systems is based on the objective of cam size, and 
is made possible for reciprocating and oscillating roller follower types. 
Optimization is not necessary for flat faced followers. 

The Zoutendijk's method of nonlinear optimization is used for the design of 
the follower parameters, and the Newton-Raphson's algorithm for the synthesis and 
analysis of linkages. 

The assembly includes the design of the assembly decemtent follower 
parameters. The type of assembler used has been given no consideration, are! the 
user needs to decide the nature of the joints in the linkage based on the type of 
the assembler he chooses to use. 

12 . Data struettres. 

The important data structures used in the system coiwjrise those for the 
representation of the trajectory, liNcages and cam-follower mechanisms. 

The trajrctory is stored as a linked list of motion sections. The struettr^ 

* 

representing each motion section contains information atout the type of space 
curve followed by the end effector in that section, ( i.e. a straight line, circular 
path, etcJ the time variation of this space curve, ( linear, harmonic, c^loidal, etc. 
); the end points of the trajectory in terms of the x and y coex^nates, and the 
starting and the terminating time instant of that section. Obviously, the 
terminating instant of the last motion section of the trajectory represents the 
cycle time needed to trace this trajectory. 

The kinematics of the end effector at small time intervals throughout the 

motion cycle, is stored in terms of the resolved components of motion in x and y 
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direction, each in an array of structures, with each structure representing the 
kinematics at one instant of time. This may also be directlu filled by the user 
without forming the linked list of the motion sections. 

The linfege is stored as a structure consisting of tte following 
information: the number of the joints with the nature of each < revolute, prismatic, 
immovable ); the flags deciding whether the vectors are of constant or variable 
lengths, and whether the angles between two vectors are fixed or variable; and a 
configuration of this linka^ for reference, to be taken as the initial point during 
analysis. These fields describe a lirtcage completely. In addition, the variables of 
analysis, and the maximum and minimun displacements of the output link are also 
stored. 

During the analysis of the linkage motion, ( that needed for cam profile 
synthesis, ) a reverse liNcage of exactly the same data type as above is formed, 
and the kinematics at small lime intervals around the time cycle is stored in an 
array of structures representing instantaneous motion as described above. The 
reason to use arrays here is to permit a quick access to the kinematics at any 
point since this will be very frequently used, especialy during the optimization of 
the follower parameters. Also, the number of divisions in i^ch the cam sutrface is 
divided is usually in the range of iOO - 120. 

The structure representing a cam mechanism consists of two fields. The 
follower parameters form one field, and are stored in arrays with predefined 
positions for each paraneter. The other field is the cam profile, ar«i this is 
stored as an array of points expressed in their x-y coordinates at small intervals 
along the cam profile. 

A table illustrating these data structures and their fields, is presented in 
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Trajectory: doubly linked list of motion sections. 

Motion section: record, 
fields: 

— motion section rsjmber (icD 

— space curve id. 

— time curve id. 

— space function parameters < array ) 

— time function parameters ( array ) 

— pointer to the previous motion section. 

— pointer to the next mdtbn section. 

Linkage: record. 

fields: 

— number of vectors in the loop. 

— joint types ( array ) 

— flag for variable link lengths. ( array ) 

— flag for variable included angles. ( array ) 

— configuration for reference. 

— information regarding the vanities of analysis. 

— maximum and minimum output displacements, 
linkage configuration: record. 

fields: 

— link lengths. < array ) 

— orientations. ( array ) 


Cam mechanism: record, 
fields: 

- follower parameters. ( array ) 

- cam profile ( array of coordinate points ) 


End effector kineealics.- docble arrau of insUntanecus kinematics 
< number of dimensions X number of divisions ) 
FcOlower kinematics: double arras of instanUneous kinematics. 

( nunkier of dimensions X number of divisions ) 


Instantaneous kinematics: record, 
fields: 

_ parameten of motion ( time / cam rotation angle ) 

— displacement. 

_ velocity. 

- acceler ation. 

r^e 7.1. Primary 


data strucUx 
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Fig. 7.1. 


73 Programme structure. 

The programme is interactive in nature. All the modules are organized in 
such a way that they form infinite lo^s embedded one inside another which can be 
exited at the will of the user. 

The programme is divided into four major nodules as follows: 

- trajectory specification. 

— synthesis of the mechanisms. 

- analysis of the mechanisms. 

— assembly. 

Each of these is described in the following sections. 

Trajectory specification. 

This module performs the following three ftHTctions. 

1. Building up of the trajectory, section by section. 

2. Determination of the resolved convonents of the end effector motion. 
SDisplay of the generated trajectory, and the end effector kinematics. 
Building tp of one motion section consists essentially of the 

identification of the space curvB, its time variation, and determining the extent of 
motion with respect to both of these. Straight lines and circular arcs are 
possible, while the time variation can te linear, harmonic, or cycloidal. The user 
can define any other path or time variation if necessary. 

Determination of the end effector kinematics is based on the curves, and 
time variations defined earlier, or direct generation of the values. 
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Synthesis. 


Synthesis comprises of two stto-modules, viz. Lirdcage synthesis and cam 
design. 

Linkage synthesis consists of generting the data for the synthesis^ and the 
actiOjal synthesis procedure. Two different modules are used to perform these two 
functions. The synthesis information collecting module gathers from the user, 
information regarding the number of variables to be synthesised, the ccKirdinated 
input and output motion values, etc. The synthesis module processes this data to 
produce the values of the linkage dimensions. 

The cam design also, comprises of two different modules. Qie generating the 
values of the follower parameters, and the other designing tt« cam profile based 
on this. In the profile design, the follower type is identified, ard a rajnber of 
points along the periphery of the cam are generated. 

There are various modules in the follower paranreter determination which 
bring about the optimization of the cam size. An input mcdule takes the data 
needed for the optimization from the user. ( tte initial values of the variables, the 
step sizes, etc., ) The actual optimization module en^loys the Zoutendilc’s method, 
and invokes two subprogranBues in this process. One is for the determination of 
the optimal search direction using the simplex algorithm, and the other is for the 
step size determination based on the direct root method. There is another module 
which continuously interacts between the Zountendijk's procedue and the follower 
optimization data, to determine the objective f motion, constraints, etc. . 
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Analysis. 


Here also, we have a subdivision in terms of the linkages, and cams. A 
linkage configuration analysis algoritttfn based on tte Newton-Raphson's technique 
is employed to determine the kinematics of the linkage at different time 
instances. Kinematics of the cam-follower system are based on a similar approach 
as that of the linkage analysis. The results of the analysis can be viewed in the 
form of various graphs, and animated displays of the lirdcage. A graph sketching 
module using the calls to the primitive graphic rescxTces in the workstation, is 
employed for this purpose. 

Assembly 

Assembly of the mechanisms is governed by two nradules. One modude 
determines the assembly dependent follower parameter, and the other determines 
the location of the cam centre. 

A chart illustrating all these modules appears in Fig. 12 . 

7.4 User interaction. 

The system is completely user interactive in nature. At every stage, the 
user has a large control in deciding the flow of the programme. Various menus 
appearing at the left hand side of the input pad, proo^it the user to take an 
appropriate action. 

I_0 

All input is through the keyboard. However, it is possible for the user to 
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TTIAJ ( trajectory specification ) 

— Editing a trajectory. <optional> 

- Motion section addition. 

— t^tion section deletion. 

— Refilling a motion section. 

— Determination of end effector kinematics. 

— Display of trajectory. 

— Path followed in space. 

— Variation of the end effector kinematics over the cycle. 


SYNTHESIS. 

- LINKAGE. 

— Specification of the nature of joints. 

— Determination of the dimensions. 

— User specified dimensions. 

— Synthesised dimensicwTS. 

— Specification of input— output coordination. 
— Solution of 2m loop closure equations. 

- CAM 

— Determination of the follower parameters. 

- User specified dimensions. 

— Optimized dimensions. 

— Specification of bwjnds, initial values etc. 
- Application of the Zoutendi jit's procedure. 
— Determination of the cam profile coordinates. 

ANALYSIS. 

-LINKAGE. 

- plots of the variation of any parameter 

of the linkage around the cycle. (Disp, vel., accln) 

- Cyclic recursive display, (superimposed) 

— Animated display. 

- CAM 

- profile 

- follower kinematics. ( disp., velocity, accln. ) 

— Pressure angle. 

- Radius of curvature. 


ASSEMBLY. 

— Location of the cam centre. 

- Determination of assembly dependent parameters, 
rw hli- 7 y Procg-aintne organization' | 



enter data from the consol, or a file, or both, intermittently. This fascilitates the 
inputting of numeric data through files, and governing the direction of flow 
through the consol. The output device is the terminal, and the graphs and other 
results are displayed on it. It is possible to direct this outjxit to a file, and a 
printer plot of the displayed figure can also be taken. 

Utmties. 

The linkage, trajectory, cam mechanism, or the kinematics of the follower and 
the end effector, can be stored and saved, so that they can be retrieved later on 
for a next run of the programme. 

Editing fascilities have teen provided for the trajectory, so that the user 
can interactively build a trajectory of the desired shape from the coded functions, 
and other user defined functions. The output is in the forms of sketctes and 
graphs, and numeric tables. All the figures can be scaled, or moved by moving the 
window. The extent of the graph may also be controlled. 

Help regarding the design procedu.'re has been inccs'porated at different 
levels of operation, and is displayed in a se;:^ate window appearing at the top of 
the screen. Helpful messages and warning are also displayed time to tiire during 
the run. 

Error handling. 

Error handling has been done to ensiF^ a correct input from the \jssr. For 
this purpose also, the programme is so organized, that it runs in an infinite loop, 
prompting the user regarding the various possible alternatives, and does not exit 

the loop unless the proper input is given. This includes taking the input amongst 

113 



the valid choices only, and to ensure a proper flow of the programme at certain 
places, for eg. the trajectory cannot be viewed unless it is built first, and the 
like. At some places, it is also possible to abort the procedtre, arwl return the 
control to a higher level. 

75 Programme control floM. 

Though the flow of the programme is governed by the user, the processing 
done from the specification of the trajectory, to the final assend^ly, is briefly 
summarized in the present section. This is also illustrated in Fig. 7.3. 

The top-most level decides the course to be taken initially. Various 
possiblities are the specification, synthesis,, analysis, or assent>ly. 

Specification of the trajectory is done in the module TRAJ. The user fills 
the linked list of the motion sections for as many motion sections as desired. 
Fascilities like adding, refilling, and deleting a motion section can te made use of 
for this. When the trajectory so formed is verified by the user, the end effector 
kinematics are calculated. The building of the trajectory is optional, and the end 
effector kinematics can be directly specified. At this stage the user has to decide 
upon the number of divisions the trajetory has to be divided into for the 
processing to be done in the succeeding modules. The end effector kinematics, and 
the trajectory can be viewed, and any modifications if necessary, can be done by 
going back to the trajectory editCK'. 

Two possible options in the s«^thesis are that of lirtcages and cams. 
Linkages are designed in the LINKAGE^SYHTHESIS module based on the erd 
effector kinematics. The extremities of the displacement are determined, and are 
displayed to fascilitate the user to specify the coordinated input and output 
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Fig. 7.3. Programme control flow. 

















values. The type of linkage depending upon the joint configtff'ation is first taken in 
form the user. This is followed by the specification of the input output 
coordinated values. The linkage is synthesised cmi the basis of these data, and the 
results are presented. If the values are not satisfactory, or a solution is not 
Possible with the specified input, the lirdcage can be resynthesised with modified 
input values. 

Cam design in the C,^M_SyNTH£SIS module comprises of the optimization of 
the follower parameters, and the profile synthesis. Direct specification of the 
follower parameters is also possible. In the optimization, the input is first 
accepted regarding the constraints, bounds on variables etc. and the module for 
the ZoutendijVs method is invoked for the optimization. The cam profile is then 
designed based on the designed follower parameters, arKi the follower kinematics. 

Various parameters of motion are then analysed in the ,4NA.ySIS module, to 
test the performance of the synthesised mechenism. Thus consists of the display of 
various plots like the velocity, acceleration of the different components of the 
meohenism, pressure angle and radius of curvatue of the cam profile, and the 
profile itself. The lirdcages can he displayed for the convilete motion cycle with 
one configur’ation superimposed on arether, or in an animated marrer. 

Any of the designed mechanisms can he resynthesised if their performance 
is not found to be satisfactory. Finally, the mechenisms in the different directions 


are assembled in the module .ASSEMBLY. 



chapter 8 . CONCLUSIONS AND 
FURTHER WORK 


8J. Ccndusions. 

In the first part of this workj issues regarding the kirematic s^thesis of a 
multi degree-of-freedran mechanism to generate a specified trajectory/ were 
studied. A design scheme for this pcsjose has been developed. 

It was seen that it is possible to generate any trajectory, differentiable, or 
non-differentiable, with these mechanisms, provided the higher derivatives of the 
components of the motion are differentiable. Cusps appear on the cam profile 
when a single CML is used to generate a non-differentiable trajectory. 

As a part of the design of the CMLs, a generalised formulation of a four 
bar planar linkage for kirematic si^thesis and analysis was done. The lir^ages 
t«are seen to serve purposes of motion magnification and type conversion < angular 
to linear arwi vice-versa ) . It was further observed, that in the absence of the 
linkages, the cam size in many cases turns out to be impracticably large. 

The major hurdle in employing thiis class of mechanisms for more than two 
degrees-of-freedom was found to be the limitation on the motion assemblers to 
maintain the independence of motion in each direction. ( i.e. for each degree of 
freedom. ) 

Motion coordination was found to play an in^jortant role in cteterminif^ the 
final trajectory generated by these mechanisms, and it was also observed that 

some flexibility in view of having modified trajectories could be_ extracted by 
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manipulating the phase angle between two cams. 

In the second part, a methcxi for the optimal design of cam-follower systems 
was proposed. It was seen that it is possible to optimize tte size of a cam subject 
to the pressure angle and radius of curvature constraints. The procedure 
developed was shown to be applicable to the reciprocating and oscillating roller 
followers. It was observed that this sort of optimization is not necessary for flat 
faced followers. 

Finally, in the third part, an interactive software was deu^loped 
incorporating the complete design of CMLs. 

82 Scope for further work. 

The present work on the CML type mechanical robots is based on the 
kinematics of the end effector and other components of the mechanism. However, 
the range of output motion itself depends on the dynamic characteristics of tte 
various moving members, particularly the assemblers. This fact can be explored by 
modeling the system desisrted, from dynamic considerations, and modifying or 
resynthesising various parameters if necessary. 

It was mentioned that the irput* output motiwi relationships in the linkages 
may help in simplifying the cam profile. In the present work, linkage design is 
based only on the input -output motion coordination. If it could be done with a view 
of simplifying the cam profile, the mechanisms can be made afWJlicable to a large 
number of seemingly difficult trajectories also. 

Multiple loop chains and 3-d cams too, could have been used in these 
mechanisms. The advantages, and limitations of using these can be studied. 

that of developing different types of motion 
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An important need is 



assemblers so that the principle suggested becomes widely applicable, especially 
to three and more degrees-of -freedom. 

In the optimization of the follower paraneters, only disc cam mechanisms 
have been considered. This principle may be applied to 3-D cams as well. The 
effects of dynamic limitations along with the kinematic constraints can also te 
included. Objective functions other than the cam size, for eg. based on the rigidity 
or dynamic performance charactehistics can be considered. 
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appendix A. THE NEWTON-RAPHSGN METHOD. 


The Newton-Raphson's method is used for the solution of nonlinear 
simultaneous equations in the linkage stffithesis and analysis. A few points are 
discussed here regarding its application. These are: 

1. Brief description of the principle. 

2. An initial point needed for the commencement of the algorithm. 

3. Criterion which decides the termination of the process. 

4. Error checking. 


Principle. 

Let 


e^ ss } . 1 Xft) — Oj 


Bf, s fft(Xi,. , . , Xjt) = 0, <A.D 

be the n equations in n variables. Let , . , x^q) be the approximate set of 
the values of the parameters, and let (AXi,. , . , Axn) te the corrections applied to 
this approximate set. Then, the equations with the new set of values become: 

fi(Xi + Axi,. , . , x„ + Ax«) = 0, 


fft(Xi+ AXj,. , . , Xn + AXff) = 0. 


(A 2) 


Expanding these by the Taylor's theorem, and neglecting the terms in his^Ter powers 
of the corrections. 


+ AXi 


'Sf 1 


+ AXb 


0 , 







* 0 . 


(A.3) 


+ AXft 


where all the coefficients of the corrections, and the constants, i^. (fi, . . . , f*) 
are ev/aluated at the current approxifnate values. These form n linear equations in 
n unknowns, and can be evaluated for the values of (Ax^,. , . , axr) by simple matrix 
transformations. . Then, the second approximation can be fouTd, and the process 
repeated till a satisfactory solution is reached. 

Initial point. 

Most of the time an initial point can be specified in the vicinty of i^ch the 
solution is expected to lie, on the basis of a^e of the following. 

1. A scaled drawing of the possible solution. 

2. Approximate calculations neglecting a few terms. 

3. Judgement based on the past exfsrience. 

The solution obtained depends on the initial point specified, and will be the 
closest possible solution from this point. The possibility of multiple solutions has 
to be explored through different sets of initial points. 

Error checking. 

This is necessary in view of the following. 




1. Divergent solutions: At times due to the nonlinearity of the equations, the 
next point obtained may lie farther off from the solution than the cirT'ent 
point, a check on the maxirmjm number of iterations is necessary to govern 
such a situation. 

2. Singularities: Many a times, the determinant of the coefficients of the 
variables in the linearized approximation may become zero. This indicates the 
existence of a singularity and a new initial point may have to be tried. 

Termination criterion. 

The solution is identified when the residual values of the equations being 
solved become nearly zero. This forms the termination criterion Ladder normal 
circumstances. The errors mentioned above may occur sometimes, and the procedure 
" need^ to be terminated under such circumstarces also. It may be retried with a 


different initial point. 


appendix B. EXPRESSIONS FOR CAM-FOLLOWER DESIGN 


The analytical expressions to map the follower motion on th« cam profile 
( as shown in Sect. 2.42 ) are: 


where 


[ ] , {R)j^ = [ M/ ] # tR)if- 

CRljc - [ Mf ] [ M, 1 , iRiif. 


(B.i> 

(B2i 


{R3lf = 

[ M/ ] , *= 

tR3lo - 
[ 3 j = 


positicMi vector of the contact point in frame Lf, 
transformation matrix from frame Lf to Cg. 
position vector of the contact point in franit- Cj 
transformation matrix from frame C to Cg. 


Analytical expressions for all the terms in the matrix I Me ] # are determined in 
section 2.42. The matrix t M/ 3 , and the vetor {R32f depend on the type of the 
follower being used, and are determined here for four different types of 
followers. Since it is possible to get analytical expressions for the x and y 
coordinates of the point on the cam profile, these are determined for eat^ of the 
case with the help of equation B2. 

Expressions for the pressure angle and radius of curvature for each case 
based on the following definitions have also been determined. 

The presstre angle 0 is defined as the angle between the direction of the 
applied force on the follower < i.e. the common normal to the cam and follower 
surfaces, e ), and the direction of the velocity vector (v) of the follower at the 



point of contact. 


Thus, 


cos ^ = g * 


(B3) 


The radius of ctrvature is defined as the reciprocal of the principal 


curvature of the cam profile at the point of contact. The CLo^ature for two 
dimensional curves is given as: [2.1 


X 


bl 

G 


<B.4) 


where N and G are the first and second ouac^atic forms for planar curves, and 
for a cam profile are given as: 


36^' see 


and, 


<B.5) 


3R aR 

° ' 36 ,' 89 , ■ 

where g is the position vector of the contact point. Both 
expressed in the same frame of reference. 


(B.6) 

and JR need to be 


Translating flat faced -follower. 

Ref. Fig. B.l. In the frame of reference L/, the velocity of the follower is 
given as: 

- CO, (i)c*Si', 0 3. (B.7) 

where, Si = the follower displacement, 

dsi 

and s.' = 

The relative velocity between the cam and follower surfaces is given as: 

Vr = Uf C (Si + So), (Si' -S - a), 0 IT (B.8a> 

The common normal to tfe cam and follower surfaces is given as: 

e = C 0, i, 0 IT. (BBb) 


Then from equation B.3, 





5 = s/ - a. 

The radius vector on the follower surface can then be written as: 

{R}2f = ( S, Si, a i IT. 

The matrix [ M/ ] ( is given as; 


[ M/ 3 


» = 


Then, from equation B2, 




COSd{ 
— sindf 
□ 

0 


1 

0 

0 

0 


sinQc 

cosdf 

0 

0 


0 

1 

0 

0 


0 

0 

1 

0 


a 

So 

0 

*1 


0 a cosdc +sosindf 

0 —a sin©« + SoCosS, 

1 0 

0 1 


= C xp, yp, Zp, 1 }, 

where, Xp = ia + S) cosQc + (so + Si) sinOc, 

Ljp = — <a + S) sin©c + <So + Si) cos©e, 

and Zp = 0. 

From equations B.3, and B.8, we get, 

= 0 . 

Also, expressing ^ in the frame C, and letting R = {R} 2 f, from equations 
and B.IO, we get, 

4 

X 

and 


(Si“ + Si + Sj)) 


(B.9) 

(B.10) 


(B.iD 


• CR>2f 

<B.i2) 


(B.i3) 

(B.i4) 
B.4, BB, 

(B.15) 


J \*6 



a 


<B.i6) 


i 


X‘ 


Translating roller follows. 


procedure is exactly similar to the ors described above^ and 
expressions involved are given here in reference to Fig. B2. 

The velocity of the follower in the frame of reference L/: 

V = CO, WfS/, 0 3. 

Si = the follower displacement, 

“ de,' 

Relative velocity between the cam and follower surfaces: 

^ = (i)c C — (Si + So + rrsin©/- ), <—5/ + rrCOsSr + all 0 3 

Common normal to the cam and follower surfaces: 

£ = £ cosOf, sin©/-, 0 3 ^. 

Condition of contact. 


—(So + Sitoos©/- + (a - SiOsinSf = 0. 

^ . -<So + Si) ^ 

_ tan < _(a _ 5 ,>) >• 

At Si = 0, ©/ lies in the third quadrant. 


Profile. 


{Rlj^ = Cr/-cos©r, (Si + rrsin©f), 0, 1 3^. 


[ M/ ] j 


1 0 0 a 
0 1 0 So 
0 0 10 
0 0 0 1 


only the 

(B.IT) 


(B.18) 

(B193 

(B20) 

(BiD 

CB22) 

<B23) 


^'7 



■ -1* 



cos6t 

-sin©« 


0 

0 


sinOf 

cos6f 

0 

0 


0 a cosQe +Sosinec 

0 —a sin9f + Socos©, 

1 0 

0 i 


Pressure angle. 


Curvature. 


= ( Xp, Wp, Zf, i }, 

Xf = rr cos(©f — ©c) + a cosSe + (sq + s^) sin©f, 

yp = Pr sin(©r — 6f) — a sin©f + (so + Sj) cos©f, 

Zp = 0. 


cin“^||sin9f 


N 

G ' 


N = — <i — ©/-Of Prii — ©/-O + (a — Si^cos©/- + (So + 

G * rr^i “ ©j'^1^ + (S{) + + (a ~ 

+ 2 rr<i — ©/■') Ka — Si')cos©r + (Sjj + Si)sin6r3. 

^ , _ Si'(a - SiO + Si"(so + si) 

= g: 5"® — 

(Sq + S^) (a — Sj^^) 

Radius of curvature. 
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Oscillating flat faced follower. 

The velocity of the follower in the frame of reference L,: (Fig. B.3> 

V = C — sine^, cos©/, 0 JT. 

6/ = the follower displacement, 


6/' 


" dec' 


Relative velocity between the cam and follower surfaces: 

^ = Ut(i + ©/O C — Ssin©/, (5 cos©/ — a / (1 -t 


Common normal to the cam and follower surfaces: 

e = ( sin©/, —cos©/, 0 }T. 

Condition of contact. 

S(1 + ©/O = a cos©/. 

Profile. 

(Rljf = C S cos©/, S sin©/, 0, i 


[ M/ 3 


t = 


—cos©/ sin©/ 
sin©/ —cos©/ 


0 

0 


(R}ic= 


0 

0 


— cos(6c+ ©/) sin(©c+ ©/) 
5in(©f+ ©/) cos(©{+ ©/) 

0 0 

0 0 


0 

0 

1 

□ 


0 

0 

-1 

0 


a 

0 

0 

1 

a COS©e 
—a sin©/ 
0 
1 


C Xf, yp, zp, 1 }, 


©/O), 0 3T 
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- cos(ef+ e^) 5 + a cosOc, 


Pressure angle. 


Curvature. 


Xf = 

yp = 

Zp = 

X = 


sin<&f+ 6/) « - a sinec, 


0 . 


Radius of curvature. 


0. 


<i + 

e/' 5 + (1 + 26, Oa sine, ' 


1 

X‘ 


Oscillating roller folloMer. 

Velocity of the follower in the frame of reference Lft 

— I 1/sin©/ + rrsin(©r + 6/) 3 
V = E 1/COsO/ + rfCos(©r + 3 

0 


s the follower displacement, 

dS, 

= 5^:- 

Relative velocity between the cam and follower surfaces: 

(1 + ©/O [ 1/sin©/ + rrsin<©r + ©/) 3 
V o (1 + 6/^) E 1/cos©^ + r^cost©/’ + ©/) 3 — a 

0 


Common normal to the cam and follower surfaces: 

e = C cos<©r + 6/), sin(6r + ©/), 0 3 ^. 
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(BAD 


(BA2) 


(BAS) 


(B.44> 


(B.45) 
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Condition of contact. 


Profile. 


Wlc= 


Pressure angle. 


Curvature. 


—(a sin6^)cos6r + ( (1 + 1/ — a cos6/ )sin6/’ = 0. 


= tan 


-1 


(a sinO/) 


( (1 + 1/ — a oos9/) 


(B.47) 
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{R}2f = f 0/ + rrCos$jr), rrsinOr, 0, 1 }T. 


[ 3 


9 = 


— cosQ^ sin©/ 

sin©/ -cos©/ 
0 0 

0 0 


-cos(©f+ ©/) sin<©<+ ©/) 
sin<©<+ ©/) cos<©f+ ©/> 


0 

0 


Xp = 

yp = 

Zp = 


0 

0 


0 

0 

1 

0 


0 

0 

-1 

0 


a 

0 

0 

i 

a cos©c 
—a sin©c 
0 
1 


<B.49) 
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^If 


<B.SD 


C Xp, yp, Zp, i }, 

— l/cos(ec+ ©/) — rp cos(©/- + ©i + ©/) + a cos©c 
1/ sin(©c+ ©/) + rr sin(6r + ©c + ©/) — a sin©c 
0. (B52) 


cos 9 = ^ s inQr 

4 I/* + Tf* + 2rr 1 /Cos©/ 

-( 1 + ^ > 

X = ffiS4) 
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<B55) 


ki « Iji + l/CosSr - a cosOf + 6,) ^ 

^2 = ^r\^f" 1/Sin6r — 6^'a cos(©^ + 

Ql + 6/^0 Q^ooE©r + Tr) — acos<§^ +©/3[] 

<BJ56> 


Radius of curvature. 
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X' 
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appendix C. EDCPREBSIQNS FGR THE ASSEMBLY OF THE CHLS. 


Expressions which are used ■for the location of the cam centre after 
determination of the follower parameters of each CH_, are derived here. Also, a 
few points as regards the assembly dependent follower parameters are mentior^d. 

C .1 Location of the cam centre. 

Location of the cam centre depends upon the combination of the follower 
type that is being used, and expressions for each are determined. Different 
combinations of oscillating and reciprocating followers are considered. For each 
case the determination forms a simple geometric point location problem. 

Two oscillating foUowers. 

The cam centre in this case is decided by the follower pivot to cam centre 
distance. Locus of the cam centre with respect to the follower pivot of each CM_, 
is a circular arc with centre at the follower pivot, and radius equal to the pivot- 
to-cam-centre distance. The cam centre should lie at the intersection of the two 
circular loci. 

If (Xi, yi) and (X2, y2) Sr's the global coordinates of the follower pivots of 
the two CMLs, and ai and a.2 the corresponding distances from the cam centre, then 
from simple geometry, the coordinates of the centre can fce given as the point of 
intersection of the two circles: 

<x — Xj)® + <y — yi)® . = ai®, tC.D 


J^3 



and 


(x - X 2 )® + <y - yj)® 


(C^ 




Ckie oscillating foUoMer and one translating foUoMer. 

The pivot-to-cam-cenlre distance of the oscillating follower, and thw offset 
of the translating follower are the parameters that decide the location of the cam 
centre. The loci of the cam centre are circular and that of a straight line fcr the 
oscillating, and the translating follower respectively. 

If (Xi, yp are the global coordinates of the follower pivot of the oscillating 
follower, and (xa, y^) that of a fixed point in space for the with a translating 
follower, and ai and aa the pivot-to-cam-centre-distance, and the offset 
repsctively, then as before, the coordinates of the centre can be given as the 
point of intersection of the circle 

I <x — Xi)* + <y — yi)* « <C5) 

and the straight line 

<y — Wz — aacos©) = tan© <x — Xa + aasin©), (C.4) 

where © is the orientation of the line of reciprocation of the follower with 
respect to the global x axis. 


Two translating followers. 

In this case, the locus of the cam centre with respect to both the followers 
is a straight line and the centre should lie at the point of intersection of these. 
The two straight lines are given as: 

<y — yi — aiCOS©i) = tan©! (x — Xi + aiSinSi), (C.5) 

<y — ya — aaCosQa) = tan©2 <x — Xa + aasinOa). <C.6) 


and 





In all these cases, the coordinates of the cam centre can te found by 
solvins the two simultaneous equations. The ntanber of solutions obtained can be 
either two, one or none. If two solutions are obtained, the one forming a more 
compact assembly may be selected. If no consistent solution is obtained with the 
two simultaneous equations, the fixed points of both the CM_s may be shifted to 
other possible locations, and the assembly tried again. It might also become 
necessary to modify the values of the offset and / or the pivot-to-cam-centre 
distance. 

CJZ Assen^ly dependent parameters. 

Since the location of the cam centre is determined after the design of the 
CMLs, it becomes necessary to transfer the temporary cam centre to that position. 
Since this might change the values of the follower displacement, two parameters, 
viz, the follower rod length, and the tte follower arm angle, are designed at a 
later stage. The values of these are determined form tl^ coordinates of the 
reference point on the follower, and the cam centre. (Fig. 2.7, 2.8). 
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apj>endix D. Some more cxjinputaiional aspects of wtalysis. 


Discussion. 

A common problem in the analysis of linkages is that one (Alains geometric 
inversions of the linkage instead of the real solulitjn at the specified value of 
the input. This happens since the lcx)p cjlosure eciuations are satisfied with that 
configuration also. Actually this inverted configuration is generally not attainable 
in the motion cycle tnless the mechanism is disassembled and reassembled in an 
inverted manner. 

However, with the Newton-RafAsco's method, we always get that solution whicA 
is in the vicinity of the initial point. Thus, if at every stage a configuration near 
to the actual solution is used as the initial point, the required solution only will 
be reached, and not its geometric inversion. For the cyclic analysis, this has the 
effect of moving the linkage through small discrete steps arocnd the entire motion 
range of the input as shown in Fig.s 3.4, 3.5. 

However, there are problems when both tte geometric inversions are nearly 
equally clcjse to the initial point. This hawjens in the case of rocker mechanisn®, 
at points near the dead center of the ir«jut link. In such mechanisms, both the 
inversions are possible with the same linkage assembly - one while approaching the 
dead center, and the other while reproaching. Physically, the path taken lp by the 
coupler and the output link are determined by their inertias. Generally, they tend 
to move in a direction in which they were previously moving. Since inertias 
cannot be modeled in the Newton-Raphson’s method, the results obtained may not be 
true to the real. There is a small compuUtional trick however, which can be used 
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to overcorne this problem. 


The fact that both, the coupler, and the output link continue to move in 
same direction as before, can be made use of. While specifying the initial point in 
the vicinity of the extremity, the values of the variables determining the notions 
of these links may be incremented in the direction of their movem«it. i.e. if Spr*// 
Scun 3^™! ®next previous, current, and next values of the variables, then 

we know that (scur — sprew) and — Scur^ have the same signs. For small 

intervals. 


or, 


scur — Sprat/ as ®next “ ®cur 
®next ® ^ Sour — ^prev- 


(3.7a) 

(3.7b) 


This can be used as the initial value for the next point, and the loop cloajre 


equations solved as before. Using this technique, we can hope to get the real 


solution. 


If however, a (*^1 occu^ at the extreme position of the input, then one 
wo(^ld not be certain about the direction in which the coupler and the output link 
would move, since Scur “ Spret/- Physically also, there wojld be no control over 
the direction of movement in such situations, since it will be determined by other 
operation dependent considerations like vibrations, clearances, etc. It is advisable 
not to use rocking mechanisms with a dwell at the input extr«nity. 
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